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PREFACE 

This  study  was  conducted  as  part  of  the  Independent  Research  Program  of  the 
Institute  for  Defense  Analyses,  under  which  significant  issues  of  general  interest  to  the 
defense  research  community  are  investigated. 
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I. 


INTRODUCTION 


RPPDM  is  a  computer  model  designed  to  solve  an  extended  version  of  the 
preallocated  preferential  defense  game  as  presented  by  Bracken,  Brooks  and  Falk  in  [1], 
The  problem  revolves  around  a  defender's  attempt  to  protect  T  targets  against  A  reentry 
vehicles.  The  targets  may  possess  different  relative  values.  (In  [1]  all  the  targets  were 
assumed  to  be  identical.)  We  will  assume  that  these  valuations  are  objective  and  that  the 
defender  is  unable  to  hide  the  values  of  particular  targets  from  the  attacker.  The  defender 
has  at  his  disposal  D  interceptors  with  which  to  destroy  the  RV's.  Each  interceptor  and  RV 
may  be  assigned  to  one  and  only  one  target.  We  assume  that  the  attacker  knows  the  total 
number  of  interceptors  but  is  unaware  of  the  specific  allocation  to  defend  each  target. 

In  [5]  and  [6]  Matheson  solved  the  case,  known  as  the  Basic  Game,  where  the 
attack  size  A  is  known  to  the  defender.  When  we  loosen  this  restriction,  however,  the 
solution  of  the  Basic  Game  no  longer  suffices,  for  the  defender  must  now  rely  on  a  single 
strategy,  which  must  be  useful  against  more  than  one  attack  size.  In  [1]  a  criterion  is 
proposed  for  a  robust  defense  which  most  closely  approximates  the  expected  outcome  of 
the  Basic  Game.  Robustnesss  problems  are  formulated  as  linear  programs  for  four  alternate 
scenarios  depicting  different  combinations  of  behavioral  assumptions  regarding  the  attacker 
and  the  defender.  The  present  model,  RPPDM,  solves  the  most  interesting  of  these 
scenarios:  the  defender  believes  that  the  attacker  can  discover  and  thus  optimize  whatever 
defense  he  chooses  to  employ,  and  the  defender  is  correct  (case  II,II  of  [1]). 

As  it  is  currently  designed,  RPPDM  may  be  used  in  either  batch  or  interactive 
mode.  In  interactive  mode,  RPPDM  prompts  the  user  for  every  piece  of  information 
necessary  to  execute  the  program.  In  batch  mode  a  separate  program  is  run  to  create  an 
input  file  for  use  with  RPPDM.  The  user  has  a  choice  of  sending  the  report  to  a  file,  the 
terminal,  or  both.  The  user  may  choose  to  solve  any  number  of  robust  defenses  for  a  given 
set  of  Basic  Game  parameters.  (This  will  be  made  clear  in  the  examples.)  The  output 
report  presents  a  summary  of  the  game  parameters,  the  optimal  strategies  for  the  Basic 
Game,  the  expected  outcome  of  the  Basic  Game,  the  robust  defense,  the  optimal  attacks 
against  the  robust  defense,  and  the  resulting  expected  target-survival  rates. 
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The  code  is  organized  around  a  main  program,  MAIN,  whose  major  function  is  to 
query  the  user  (or  read  from  a  file)  for  specific  parameters  needed  to  set  up  the  problem, 
such  as  the  number  of  interceptors,  the  number  of  targets,  etc.  MAIN  then  uses  that 
information  to  call  on  a  series  of  subroutines  to  conduct  the  appropriate  numerical 
caluclations  and  to  generate  the  solution  reports. 

The  subroutines  are  organized  into  three  categories: 

1)  PU  subroutines 

2)  LP  subroutines 

3)  REPORT  subroutines 

The  PU  subroutines  include  SIMAT1,  SIMAT2,  SEQAT1,  and  SEQAT2.  They 
are  used  to  generate  the  PU's  associated  with  particular  combinations  of  attack  and  defense 
methodologies.  Only  one  of  these  subroutines  is  used  in  any  one  run  of  RPPD. 

The  LP  subroutines  include  BG,  YROUBUST  and  XROBUST.  These  subroutines 
are  the  heart  of  the  model.  They  set  up  the  linear  programming  equivalents  of  the  problem 
in  a  format  that  can  be  accepted  by  XMP,  a  linear  programming  system.  BG  is  responsible 
for  returning  the  optimal  attacker  and  defender  strategies  the  expected  outcome  (game 
value)  of  the  Basic  Game.  YROBUST  calculates  the  robust  defense.  XROBUST  solves 
for  the  optimal  attacks  against  the  robust  defense  and  the  resulting  expected  survival  rates. 

The  REPORT  subroutines  include  SUMMARY,  STPRINT,  YPRINT,  VPRINT, 
ALPRINT,  ALYRPRINT,  AL VPRINT,  and  RVINTCOUNT.  They  are  used  to  produce 
the  output  reports. 

The  code  is  written  in  FORTRAN-77  to  run  on  a  VAX  8600  computer. 

Note  that  computer  programs  which  solve  the  Basic  Game  are  documented  in  [3], 
14],  and  [7]. 

II.  MAIN  PROGRAM 

Program  MAIN  acts  as  a  coordinator  between  the  user  and  the  subroutines  where 
the  actual  "work"  takes  place.  MAIN  first  asks  the  user  whether  interactive  or  batch  mode 
is  to  be  used.  If  batch  was  chosen  it  asks  fo.  the  name  of  the  input  file.  The  flow  of 
MAIN  may  then  be  broken  roughly  into  five  stages: 

1)  Generate  the  PU’s 

2)  Solve  the  Basic  Game 


vV  v.V v 


3)  Produce  the  solution  report  for  the  Basic  Game 

4)  Solve  the  robust  game 

5)  Produce  the  solution  report  for  the  robust  game 

During  each  stage,  MAIN  will  prompt  the  user  (or  read  from  file  in  batch  mode)  for  the 
additional  information  needed  to  process  that  stage.  We  will  discuss  each  of  these  stages 
separately,  focusing  on  the  needed  inputs  and  the  subroutines  called. 

Stage  1;  Generate  the  PIJ's 

RPPDM  is  equipped  to  handle  five  combinations  of  attack  and  defense 
methodologies: 

A)  Simultaneous  Attack  with  One  Shot 

B)  Simultaneous  Attack  with  Shoot-Look-Shoot 

C)  Sequential  Attack  of  Unknown  Size  with  One  Shot 

D)  Sequential  Attack  of  Unknown  Size  with  Shoot-Look- Shoot 

E)  Sequential  Attack  of  Known  Size  with  Shoot-Look-Shoot 

Depending  upon  the  answers  to  two  questions  (three  if  the  attack  is  sequential),  MAIN  will 
call  on  one  of  four  subroutines,  SIMATT1  for  A,  SIMAT2  for  B,  SEQAT1  for  C  and  D, 
and  SEQAT2  for  E,  to  generate  the  appropriate  set  of  PIJ's.  For  each  of  the  PIJ 
subroutines  the  following  inputs  are  necssary:  the  maximum  number  of  RVs  that  may 
defend  a  single  target  (S),  the  failure  probability  of  the  RV’s  (PFA)  and  that  of  the 
interceptors  (PFD).  For  cases  involving  shoot  look  shoot  defenses,  two  failure 
probabilities  for  the  interceptors  are  necessary  (PFD1  and  PFD2  rather  than  PFD). 

Stags  2;  Solve  the  Basic  Game 

Main  prompts  for  the  minimum  attack  size  (MINRV),  the  maximum  attack  size 
(MAXRV),  and  the  attack  size  increment  (INCRV).  These  values  are  then  used  to  calculate 
the  actual  set  of  attack  sizes  (RV)  and  the  number  of  attack  sizes  (A).  The  number  of 
interceptors  (INT)  and  the  total  number  of  targets  (TARGETS)  are  also  needed.  In 
addition,  if  there  are  more  than  one  type  of  target,  the  relative  value  VTYPE)  and  the 
number  of  targets  (NTAR)  belonging  to  each  type  must  also  be  entered.  The  value  and  the 
size  of  each  target  type  is  then  calculated  as  the  fraction  of  all  the  targets.  Subroutine  BG  is 
called  to  solve  the  Basic  Game.  BG  returns  the  attracker's  minimax  strategy  in  XBG,  the 
defender's  minimax  strategies  in  YBG,  and  the  game  values  in  VBG. 


MAIN  calls  on  the  REPORT  subroutines  to  print  out  the  results  of  the  Basic  Game 
on  the  selected  output  device(s).  SUMMARY  provides  a  listing  of  the  problem’s 
parameters  as  entered.  STPRINT  prints  out  XBG  and  YBG,  VPRINT  prints  out  VBG, 
ALPRINT  prints  out  allocation  tables  based  on  XBG  and  YBG,  AL VPRINT  prints  out  the 
expected  number  of  targets  that  will  survive  for  each  target  type,  and  RVINTCOUNT 
prints  out  the  numbers  RV's  and  interceptors  assigned  to  each  target  type.  (See  Section 
VII,  "Notes  on  Output"  for  more  details.) 


Stage  4:  Solve  the  Robust  Game 

MAIN  prompts  for  the  number  of  robust  strategies  desired.  Stages  4  and  5  are 
repeated  until  all  of  them  have  been  solved 

MAIN  prompts  for  a  lower  and  an  upper  bound  on  the  attack  sizes  for  which  a 
robust  defense  is  to  be  found.  Both  bounds  must  lie  in  the  set  of  attack  sizes  defined  for 
the  Basic  Game  during  stage  2.  Subroutine  YROBUST  is  then  called  to  solve  for  the 
robust  defense  (YII).  With  the  robust  defense  defined,  MAIN  calls  XROBUST  to  find  the 
optimal  attacks  (XII)  against  YII  and  the  resulting  expected  survival  rates  (VII). 

Stage  5-  Produce  the  solution  report  for  the  Robust  Game 

MAIN  calls  on  YPRINT  to  print  out  YU,  STPRINT  to  print  out  XII,  VPRINT  to 
print  out  VTI,  AL  YPRINT  to  print  out  the  allocation  table  for  YII,  ALPRINT  to  print  out 
the  allocation  table  for  XII,  AL  VPRINT  to  print  out  the  expected  number  of  targets  that  will 
survive  for  each  target  type,  and  RVINTCOUNT  to  print  out  the  numbers  of  RV's  and 
interceptors  assigned  to  each  target  type. 

III.  PU  SUBROUTINES 

The  PU  subroutines  SIMAT,  SIMAT2,  SEQAT1,  and  SEQAT2  are  called  on  to 
generate  the  appropriate  set  of  PIJ’s  associated  with  particular  attack  and  defense 
methodologies. 

SIMAT1  is  designed  to  generate  the  PU's  when  the  attack  is  simultaneous  and  the 
defense  has  only  one  opportunity  to  intercept  the  RV’s. 


SIMAT2  applies  when  the  attack  is  still  simultaneous,  but  the  defense  has  two 
opportunities  to  intercept  the  RV. 

SEQAT1  applies  when  the  attack  is  sequential  of  unknown  size  and  the  defender 
has  one  or  two  chances  to  intercept  the  RV's.  The  first  case  is  accomplished  by  passing 
1.0  for  the  failure  rate  of  the  first  interceptor  salve  (PFD1)  and  assigning  PFD  to  PFD2.  In 
both  cases  the  defender  does  not  know  the  number  of  RV's  that  will  arrive  at  a  particular 
target. 

SEQAT2  applies  when  the  attack  is  sequential,  the  defender  has  shoot-look-shoot, 
and  the  defender  knows  (after  the  attack  begins)  how  many  RV’s  will  engage  him  at  a 
given  target. 

The  resulting  PIJ's  are  returned  to  MAIN  in  the  two  dimensional  array  P.  The  first 
index  of  P  is  the  number  of  RV's  that  attacks  a  target,  and  the  second  index  is  the  number 
of  interceptors  that  defends  a  target.  Thus  P(I,J)  is  the  probability  that  a  target  will  survive 
when  attacked  by  I  RV's  and  defended  by  J  interceptors.  See  Appendix  H  for  discussions 
of  the  concepts  on  which  these  PU  subroutines  are  based. 

IV.  LP  SUBROUTINES 

The  purpose  of  BG,  YROBUST,  and  XROBUST  is  to  set  up  the  appropriate  linear 
programs  for  processing  by  XMP.  In  the  following  subsections  we  will  examine  each 
individually  in  terms  of  the  linear  programs  they  solve.  Please  consult  program  listing  in 
the  Appendix  and  XMP  documentation  for  details  regarding  the  mechanics  of  the 
interactions  between  BG,  etc.,  and  XMP.  In  the  linear  programs  that  follow,  all  variables 
are,  or  can  be  assumed  to  be,  non-negative. 

m 

BG  is  designed  to  solve  the  Basic  Game.  A  linear  program  of  the  Basic  Game  may 
be  set  up  as  follows  (see  Appendix  G  for  the  derivation): 
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ktl  J 

subject  to 
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The  problem  is  solved  first  for  RV(1),  the  smallest  attack  size,  and  the  solutions 
stored  in  VBG,  XBG  (the  dual  variables  of  the  inequality  constraints)  and  YBG.  Then 
RV(1)  is  substituted  by  RV(2)  and  the  problem  reoptimized.  BG  keeps  substituting  RV's 
until  all  of  the  attack  sizes  are  solved. 

XBG  and  YBG  are  passed  back  to  MAIN  as  three  dimensional  arrays.  The  first 
index  designates  the  attack  size,  the  second  index  is  the  target  type,  and  the  third  index  is 
the  number  of  RV's  (for  XBG)  or  interceptors  (for  YBG).  XBG(A,K,I)  and 
YBG(  A,K,J),  are,  respectively,  the  fraction  of  type  K  targets  attacked  by  I  RV's  and  the 
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fraction  of  type  K  targets  defended  by  J  interceptors  for  the  basic  game  with  .he  A'th  attack 
size  (RV(A)). 

VBG  is  passed  back  to  MAIN  as  a  one  dimensional  array  with  the  index 
designating  the  attack  size.  VBG(A)  is  the  game  value  for  the  A'th  attack  size. 


YROBUST 

YROBUST  is  designed  to  solve  for  the  robust  defense.  A  linear  program  of  the 
robust  game  may  be  set  up  as  follows  (see  Appendix  G  for  the  derivation’): 


max  P 

Y(k,j),s(A,k),t(A) 


subject  to 


VF(k)  X 
j=0 
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VFW  % 
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The  Y  that  solves  this  program  is  stored  m  YD  and  passed  back  to  MAIN  as  a  two 
dimensional  variable.  The  first  index  is  target  type,  and  the  second  index  is  the  number  of 
interceptors  assigned  to  a  target  YII(KJ)  is  the  fraction  of  type  K  targets  defended  by  J 
interceptors  under  the  robust  defense. 


XROBUST 


Given  the  robust  defense,  YII,  XROBUST  finds  the  optimal  attack,  XII,  and  the 
resulting  expected  survival  rate  using  the  following  LP: 


NTYPE  R 

VH  =  min  £  ■ 

X(k,i)  k=l  i=0 


subject  to 


£P(i,j)  YII(k,j) 
j=0 


X(k,i) 


R 


X  X(k,i)  =  l 
i=0  _ 


1  £k< NTYPE 


NTYPE  R 


I  I 


k=l  i=0 


i-  NF(k)  X(k,i) 


=  RV/TARGETS 


After  the  LP  is  solved,  and  the  X  that  yields  VII  is  stored  in  XII  for  every  attack 

size. 

XII  and  VII  are  passed  back  to  MAIN  as,  respectively,  two  and  one  dimensional 
variables.  The  indices  are  the  same  as  their  counterparts  in  the  Basic  Game,  XBG  and 
VBG. 

V.  REPORT  SUBROUTINES 

The  report  subroutines  generate  formatted  output  for  RPPDM.  The  SUMMARY 
subroutine  creates  a  copy  of  the  parameters  of  the  game  as  specified  and  send  it  to  the 
output  device  selected.  The  STPRINT  subroutine  prints  out  attacker  and  defender 
strategies  that  are  indexed  by  attack  sizes  (i.  e.,  XBG,  YBG,  and  XII).  The  YPRINT 


subroutine  prints  out  the  robust  defense.  The  VPRINT  subroutine  prints  out  the  expected 
survival  rate  for  the  chosen  range  of  attack  sizes.  ALPRINT  prints  out  target  allocation 
tables  for  attacker  and  defender  strategies  indexed  by  attack  sizes.  ALYPRINT  prints  out 
target  allocation  tables  for  the  robust  defense.  AL VPRINT  prints  out  the  expected  numbers 
of  targets  that  will  survive  for  each  target  type.  RVINTCOUNT  prints  out  the  missile  (RV 
or  interceptor)  allocation  by  target  types. 

VI.  NOTES  ON  INPUT 

A  separate  program,  BATCH,  is  provided  to  facilitate  the  use  of  batch  mode.  It  is 
an  interactive  program  that  creates  an  input  file  that  may  be  read  by  RPPDM.  It  will  query 
the  user  for  every  piece  of  information  that  is  needed  to  run  RPPDM  (in  the  same  order  as 
the  interactive  mode),  and  write  the  information  to  a  sequential  file  named  by  the  user. 
Thus  when  a  user  seeks  to  run  a  series  of  problems  with  changes  in  only  a  small  number  of 
the  parameters,  he  or  she  could  create  the  first  problem  using  BATCH,  then  use  an  editor  to 
modify  the  input  file  for  the  first  problem  to  fit  the  data  for  the  other  problems.  We 
demonstrate  this  procedure  in  the  second  example  of  Section  VUI. 

VII.  NOTES  ON  OUTPUT 

In  any  given  run  of  RPPDM  a  minimum  of  9  tables  are  printed,  and  for  each  robust 
defense  desired,  an  additional  8  tables  are  added.  (For  runs  involving  targets  of  equal 
value,  the  numbers  are  7  and  6,  since  the  RV  and  interceptor  allocation  tables  are 
eliminated.) 

The  first  table  is  a  table  of  the  basic  parameters  that  define  the  Basic  and  the  Robust 
Games. 

The  next  eight  tables  list  the  results  of  the  Basic  Game.  The  first  two  of  these  tables 
print  out,  respectively,  the  attacker  and  defender  strategies  indexed  by  target  types  and  then 
by  attack  sizes.  For  each  attack  size  A,  there  are  N  entries  which  define  the  minimax 
strategy  for  a  given  target  type  K.  The  i'th  entry  corresponds  to  the  fraction  of  type  K 
targets  that  will  be  assigned  i-1  RV's  or  interceptors.  Each  row  contains  a  maximum  of  10 
entries.  The  first  row  will  include  the  first  10  entries,  the  second  row  the  next  10  entries, 
etc.  Therefore,  an  entry  in  the  third  row  and  sixth  column  is  the  fraction  of  targets  (of  that 
target  type)  to  receive  25  RV’s  or  interceptors. 


The  third  table  of  the  Basic  Game  is  the  game  value,  or  the  expected  survival  rate 
associated  with  the  minimax  strategies  of  the  Basic  Game.  Each  entry  is  the  fraction  of  the 
"total  value"  that  is  expected  to  survive.  Thus  .3567  means  that  the  expected  value  of  all 
the  targets  that  survive  the  attack  is  35.7%  of  the  original  value  before  the  attack. 

The  fourth  and  the  fifth  tables  are  just  a  variation  on  the  strategy  tables  discussed 
earlier.  These  target  allocation  tables,  instead  of  using  fractions,  give  the  actual  numbers 
that  correspond  to  the  minimax  strategies.  The  allocation  tables  are  first  subdivided  by 
attack  sizes  and  then  by  target  types.  For  each  attack  size  the  strategy  for  all  the  different 
target  types  are  listed  together.  The  Kth  entry  now  corresponds  to  the  number  of  targets 
(of  the  target  type  under  consideration)  that  should  be  assigned  K- 1  missiles. 

The  seventh  table  lists  the  number  of  targets  of  each  target  type  that  is  expected  to 
survive  in  the  Basic  Game. 

The  eighth  and  ninth  tables  lists  the  RV  and  interceptor  allocations  by  target  type. 
Each  entry  corresponds  to  the  total  number  of  missiles  assigned  to  all  the  targets  of  that 
target  type.  (As  noted  above,  these  tables  are  not  created  when  there  is  only  one  type  of 
target.) 

The  next  eight  tables  correspond  to  the  first  Robust  Game  (if  any).  The  tables  are 
printed  out  and  the  entries  defined  in  an  identical  way  to  those  in  the  Basic  Game  except 
that  the  defender's  strategy  and  target  allocation  tables  are  not  indexed  by  attack  sizes.  The 
first  table  lists  the  robust  defender's  strategy,  the  second  the  optimal  attacker  strategies 
against  the  robust  defense,  the  third  the  expected  value  of  teh  robust  game,  the  fourth  and 
fifth  the  target  allocation  tables  for  the  defender  and  the  attacker  (respectively),  the  sixth  the 
expected  number  of  targets  surviving,  by  target  type,  and  the  seventh  and  eighth  the  missile 
allocation  tables  for  the  attacker  and  the  defender  in  this,  the  first  Robust  Game. 

If  another  robust  defense  is  called  for,  then  eight  more  tables  will  be  printed  in  the 
same  order  as  the  eight  associated  with  the  first  Robust  Game. 

VIII.  HOW  TO  USE  RPPDM 

The  following  is  a  list  of  the  FORTRAN  files  in  RPPDM  and  their  contents: 
File  Contents 


RPPDM.FOR 


MAIN 


PU1.FOR 

SIMAT1 

PU2.FOR 

SIMAT2 

PU3.FOR 

SEQAT1 

PU4.FOR 

SEQAT2 

LP1JFOR 

BG 

LP2.FOR 

YROBUST 

LP3.FOR 

XROBUST 

REPORTl.FOR 

SUMMARY, STPRINT 
YPRINT,  VPRINT 

REPORT2.FOR 

ALYPRINT 

ALPRINT,  AL VPRINT 

REOIRT3.FOR 

RVINTCOUNT 

BATCH.FOR 

BATCH 

In  addition,  the  XMP  linear  programming  package  is  necessary  for  the  use  of  RPPDM. 

Together  with  the  XMP  subroutines,  all  the  files  except  BATH.FOR  must  be 
compiled  and  then  linked  into  an  executable  image.  BATCH.FOR  should  be  compiled  and 
linked  separately. 

The  following  examples  illustrate  typical  runs  of  the  model.  In  the  first  example  the 
user  wishes  to  examine  the  robust  defense  under  simultaneous  attack  and  one  shot  defense 
with  only  one  type  of  target.  In  the  second  example  the  user  wants  to  find  two  different 
robust  defenses  under  a  more  complicated  combination  of  attack  and  defense 
methodologies  with  3  different  types  of  targets.  In  the  first  case  he  or  she  interactively 
provides  the  model  with  all  the  specifications  needed.  In  the  second  case  he  or  she  will 
create  an  input  file  with  BATCH  and  then  use  batch  mode  with  RPPDM. 

In  order  to  solve  the  Basic  Game,  RPPDM  needs  the  attack  and  defense 
methodologies,  the  RV  and  interceptor  failure  rates,  the  minimum  and  maximum  attack 
sizes,  the  attack  size  increment,  the  number  of  interceptors,  the  total  number  of  targets,  the 
relative  value  of  each  target  type,  and  the  number  of  targets  in  each  target  type.  For  the 
Robust  Game,  RPPDM  needs  to  know  the  minimum  and  maximum  attack  sizes  in  the 
range  of  attack  sizes  for  which  the  robust  defense  is  to  be  found;  both  must  be  in  the  set  of 
attack  sizes  for  the  Basic  Game.  The  same  increment  from  the  Basic  Game  is  used  here. 

These  sample  runs  are  shown  exactly  as  they  would  have  appeared  on  a  terminal 
screen  except  for  two  items:  the  '*****USER  SAYS  >'  segment  that  preceed  every  user 
entry  and  the  comments  outlined  by  slashes.  They  are  provided  to  highlight  and  clarify 
certain  situations. 
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XMAPS. . .YOU  HAVE  ROOM  FOR  2987  N0N-2ER0S  IN  THE  BASIS  FACTORS 
YOU  COULD  REDUCE  REAL  MEMORY  FROM  1000  TO  591 


SOLUTION  REPORT  FOR  THE  ROBUST  GAME 


EXPECTED  TARGET  SURVIVAL  RATE  WITH  THE  ROBUST  DEFENSE 


the  optimal  attack  allocation  against  the  robust  defense 
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TARGET  TYPE  3:  30.00X  OF  TOTAL  TARGETS,  WITH  48.84X  OF  TOTAL  VALUE 
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HE  ATTACKER ' S  BASIC  GAME  RV  ALLOCATION  BY  TARGET  TYPE 
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APPENDIX  A 
PROGRAM  MAIN 


"1 


•a 


PROGRAM  MAIN 


MAIN  will  solve  for  the  basic  game  and  the  Ca6e  II , II  as 
outlined  in  Bracken,  Brooks,  and  Falk  [1}.  Elth  batch  or 
Interactive  modes  may  be  employed.  IMP  is  the  linear  programming 
package  used. 

Parameter  statements  for  establishing  limits  on  the  problem's 
size 

PARAMETER  (MAXNAT-2S,  MAXR-30 ,  MAXS-30,  MAXNTYPE-7) 


MAXNAT  THE  MAXIMUM  ALLOWABLE  NUMBER  OF  ATTACH  SIZES 

MAXR  THE  MAXIMUM  ALLOWABLE  'R' 

MAXS  THE  MAXIMUM  ALLOWABLE  'S' 

MAXNTYPB  THE  MAXIMUM  ALLOWABLE  TYPES  OF  TARGETS 


Variable  type  declarations  for  the  basic  game 

REAL  XBGCMAXNAT,  MAXNTYPB,  O-'MAXR) 

REAL  YBG( MAXNAT,  MAXNTYPB.  0:MAXS),  VBG( MAXNAT) 

REAL  VTYPBC MAXNTYPB),  VPRAC ( MAXNTYPB ) .  NFRACC MAXNTYPB) 
REAL  P(0:MAXR,  0:MAXS).  TS ( MAXNAT ,  MAXNTYPB) 

INTEGER  R,  S . MAXRV . MINRV . INCRV , INT , TARGETS 
INTEGER  RV(MAXNAT) , A, NT AR( MAXNTYPB) 

INTEGER  TKXBG( MAXNAT. MAXNTYPB) , TMYBGC MAXNAT , MAXNTYPE) 


THE  ATTACKER ' S  OPTIMAL  BASIC  GAME  STRATEGY 
THB  DEFENDER'S  OPTIMAL  BASIC  GAME  STRATEGY 
THE  GAME  VALUE 

THE  BXPBCTBD  NUMBER  OF  TARGETS  (SURVIVING) 

THE  PIJ  MATRIX 

THB  MAXIMUM  NUMBER  OF  RV’S  AT  A  SINGLE  TARGET 
THE  MAXIMUM  NUMBER  OF  INTERCEPTORS  AT  A  SINGLE 
TARGET 

THB  MAXIMUM  ATTACK  SIZE 

THE  MINIMUM  ATTACK  SIZB 

THB  ATTACK  SIZE  INCREMENT 

THE  NUMBER  OF  INTERCBPTORS 

THB  NUMBER  OF  TARGETS 

THE  ATTACK  SIZES 

THE  HUMBER  OF  ATTACK  SIZES 

THE  RELATIVE  VALUES  OF  THE  TARGET  TYPES 

THB  NUMBER  OF  TARGETS  FOR  EACH  TYPE 

THE  FRACTION  OF  TOTAL  VALUE  FOR  EACH  TYPB 

THE  FRACTION  OF  TOTAL  TARGETS  FOR  BACH  TYPB 

TARGBT  EQUIVALENTS  OF  EXPECTED  ' VBG '  AND  'VII' 

RV  ALLOC.  FOR  BASIC  GAME  BY  TARGET  TYPES 
INTERCEPTOR  ALLOC.  FOR  BASIC  GAME  BY  TARGET  TYPES 


MAXRV 

MINRV 

INCRV 

INT 

TARGBT S 
RV 

A 

VTYPB 

NTAR 

VPRAC 

NFRAC 

TS 

TMXBG 

TMYBG 


c 

c 

o 


Variable  type  declarations  for  Case  XI. II 


REAL  YIKMAINTYPS.  0:MAXS) 

REAL  XII(MAXNAT , MAXNTYPE , 0 : MAXR) ,  VII(KAXNAT) 

REAL  VBGR(MAXHAT) 

INTEGER  MAXRRV,  MINRRV 
INTEGBR  RRV ( MAIN AT ) .  AR 

INTEGER  TMXII(MAXNAT,MAXNTYPB) ,  TMYIICMAXNAT, MAXNTYPE) 
o 

o  •»•••««• 

O  YII 

c  III 

C  VII 

C  VBGR 

C  MAXRRV 

o  MINRRV 

C  RRV 

C  AR 

C  TMXII 

O  TMYII 

c  ••••«••• 

o 
o 

INTEGER  TER . IN , PILE , 00T . OOT1 . OUT2 , QPR 
CHARACTER* 1  ANSWER .  NAME* 80,  FILBOUT*12,  FILEIN*12 


THB  ROBUST  DEFENDER'S  STRATEGY 
THE  OPTIMAL  ATTACKER'S  STRATEGY  AGAINST  YII 
THE  EXPECTED  SURVIVAL  RATE  GIVEN  YII  «JTD  XII 
THE  GAME  VALUB  FOR  ROBUST  ATTACK  SIZB  RANGE 
THB  MAXIMUM  ATTACK  SIZE  IN  THE  ROBUST  RANGE 
THB  MINIMUM  ATTACK  SIZE  IN  THE  ROBUST  RANGE 
THB  ATTACK  SIZES  IN  THE  ROBUST  RANGE 
THB  NUMBER  OF  ATTACK  SIZES  IN  THE  ROBUST  RANGE 
RV  ALLOC.  FOR  BASIC  GAMB  BY  TARGET  TYPES 
INTERCEPTOR  ALLOC.  FOR  BASIC  GAME  BY  TARGET  TYPES 


a 

c 

o 

o 

o 

o 

0 

o 

o 

o 

o 

o 

c 

o 


o 

C 

o 

o 

C 

0 

o 


TER  THB  I/O  UNIT  NUMBER  FOR  THB  TERMINAL 

IN  THE  I/O  UNIT  NUMBER  FOR  THE  INPUT  DEVICE 

FILE  THB  I/O  UNIT  NUMBER  FOR  AN  OUTPUT  FILE 

OUT  THE  I/O  UNIT  NUMBER  FOR  THE  OUTPUT  DEVICE 

QPR  THB  I/O  UNIT  NUMBER  FOR  DISPLAY  OF  PROMPTS 

ANSWER  Y'  OR  N'  ANSWER  TO  SELECTED  QUESTIONS 

NAME  HEADER  FOR  OUTPUT 

FILBOUT  NAME  OF  AN  OUTPUT  FILE 


Set  up  oommon  statement  for  input,  IOIN,  error  output,  IOERR, 
results  output,  IOLOG  (XMP  only) 

COMMON / I0_UHIT / IOIN ,  IOERR,  IOLOG 

IOIN-8 

IOERR* B 

IOLOG-B 

Set  TER  to  default  I/O  unit  number  for  terminal 


TBR-5 


Set  FILE  to  TSR+1 
PILB-TBR+1 


A-2 


oo  oooooo  ooo 


I 

I 


f 

I 


c 

C 

o 


o 

o 

c 


Set  IN  to  the  I/O  unit  number  for  the  input  file/device 

WRITE ( TER , ' ( A ) ' )  ' OTHERE  ARE  TWO  INPUT  OPTIONS:' 

WRITE (TER, ' (A) ' )  '0  1)  TERMINAL  (INTERACTIVE)' 

WRITE (TER, ' (A) ' )  '  2)  FILE  (BATCH)' 

WRITE ( TER , ' ( A ) ' )  ' OPLEASE  ENTER  THE  NUMBER  OF  THE  DESIRED  OPTI 
ION  ?' 

READ(TER, * )  NIN 
IF  (NIN  .NE.  1)  THEN 

WRITE(TER , ' (A) ' )  'OPLEASE  ENTER  THE  NAME  OF  THE  INPUT  FILE 
l(OF  LESS  THAN  12' 

WRITE(TER, ' (A) ' )  '  CHARACTERS  INCLUDING  THE  EXTENSION)  ?' 
READ (TER,  '(A)')  FILEIN 


Open  file  for  file  input 

OPEN(UNIT-TER+2 ,  FILE-FILEIN.  STATUS- ' OLD ' ) 

IN-TER+2 

Open  Junk  file  for  unnecessary  output 

OPEN(UNIT-TER+3 ,  FILE- ' JUNE . DAT ' .  STATUS- ' SCRATCH ' ) 
QPR-TER+3 

ELSE 

gPR-TER 
IN-TER 
END  IF 

Select  output  option 


write(QPR, ' (A) ' )  ' lYou  have  three  options  for  the  output  of  th 

Xo  r68Ul1)8  1 

WRITE(gPR, ' (A) ' )  '0  1)  TERMINAL  only' 

WRITE(gPR, '(A)')  '  2)  FILE  only' 

WRITE(gPR, ' (A) ' )  '  3)  TERMINAL  and  FILE' 

WRITE (QPR, ' (A) ' )  'OPlease  enter  the  number  for  the  desired  opt 
lion  ?' 

READ(IN.  FMT-* )  NOUT 


Ask  for  file  name,  if  neoessary 


IF  (NOUT  .NE.  1)  THEN 

WRITB(gPR, ' (A) ' )  'OPlease  type  in  the  desired  file  name  (o 
If  less  than  10  characters' 

WRiTB(gPR, ' (A) ' )  *  Including  the  extension)  ?' 

READ(IN, '(A)')  FILEOUT 
0P8N( UNIT-FILE,  FILE-FILBOUT ,  STATUS- ' NEW ’ ) 

ENDIF 


Input  R  and  S 

WRITE(UNIT-QPR,  FMT- ’ (A , 12 , A) ' )  ' OThe  MAXIMUM  number  of  RV’s 
l(up  to  ' , MAXR, ' )  at  a  single  target  ?  ' 

RBAD(UNIT-IN,  FMT-*)  R 


A-  3 


TRITE (UNIT-QPR,  FMT- ' (A, 12, A) ' )  ' OThe  MAXIMUM  number  of  INTERC 
1BPTORS  (up  to  '.MAXS,')  at  &  single  target  ?  ' 

RSAD(UNIT-IN,  FMT-*)  S 

Input  failure  rates  for 

Select  attach  methodology 

TRITE (UNIT-QPR,  FMT-'(A)')  ' OSelect  one  of  the  following  attac 
lk  methodologies :  ' 

TRITE (UNIT-QPR,  FMT- '(A)')  '0  1)  SIMULTANEOUS  ATTACK' 

TRITE (UNIT-QPR,  FMT- ‘ (A) ' )  '  2)  SEQUENTIAL  ATTACK ' 

TRITB ( UNIT-QPR ,  FMT-'(A)')  'OPlease  input  the  number  of  the  de 
lsired  attaok  ?  ' 

RBADCUNIT-IN,  FMT-*)  NATTYPB 

Input  the  failure  rate  for  the  RV''s 

TRITE (UNIT-QPR.  FMT- '(A)')  'OThe  FAILURE  rate  of  the  RV"s  ?  ' 
READ (UNIT -IN,  FMT-*)  PFA 

Select  defense  methodology 

trite (UNIT-QPR.  FMT- ' (A) ' )  'OSelect  one  of  the  following  defen 
lse  methodologies: 

TRITE (UNIT-QPR,  FMT-'(A)')  '0  1)  ONE  SHOT' 

TRITE (UNIT-QPR,  FMT- '(A)')  '  2)  SHOOT  LOOK  SHOOT' 

TRITE (UNIT-QPR.  FMT- ' (A) ' )  'OPlease  input  the  number  for  the  d 
leslred  option  ?  ' 

READ( UNIT-IN ,  FMT-*)  NDFTYPB 

If  the  the  attack  is  sequential,  find  out  whether  the  defender  knows, 
after  the  attaok  begins,  the  number  of  RV's  slated  for  each  target. 

IF  (NATTYPB  .EQ.  2)  THEN 

TRITB (UNIT-QPR,  FMT-'(A)')  OIs  the  defender  aware,  after 
lthe  attack  begins,  of  the  number' 

TRITE(UNIT-QPR ,  PMT-'(A)')  '  of  RV's  slated  for  each  targ 
let  (Y  or  N)  ?' 

READ(UNIT-IN,  FMT- '(A)')  ANSTER 
IF  (ANSTER  .BQ.  'Y')  THEN 
IF  (NDFTYPB  .BQ.  1)  THEN 

TRITB(UNIT-QPR ,  FMT- ' (A) ' ) ' ONOTB :  This  soenario  is  equ 
llvalent  to  one  with  a  simultaneous  attaok. ' 

NATTYPB- 1 
BNDIF 
BNDIF 
BNDIF 

Find  the  failure  rates  for  the  lnteroeptors 
IF  (NDFTYPB  .EQ.  1)  THEN 

TRITE (UNIT-QPR,  FMT-  (A)')  OThe  FAILURE  rate  of  the  inter 
loeptors  ?  ' 


RBAD( UNIT-IN,  FKT- * )  PFD 

ELSE 

WRITB(UNIT-QPR,  FMT-'(A)')  'OThe  FAILURE  rate  for  the  firs 
It  salvo  interceptors  ?  ' 

RSAD(UNIT-IN,  FKT- • )  PFD1 

WRITE (UNIT-QPR,  FMT-'(A)')  'OThe  FAILURE  rate  for  the  eeco 
lnd  salvo  interceptors  ?  ' 

READ( UNIT-IN,  FMT-* )  PFD2 
ENDIF 

Generate  the  appropriate  pij's 

IF  (NATTYPB  . EQ.  1  .AND.  NDFTYPE  . EQ.  1)  THEN 

Simultaneous  attach  without  SLS 

CALL  SIMAT1(MAER,  MAES,  PFA,  PFD,  R,  S.  P) 

BLSB  IF  (NATTYPB  .BQ.  1  .AND.  NDFTYPE  .  EQ.  2)  THEN 

Simultaneous  attach  with  SLS 

CALL  SIMAT2(MAER ,  MAES,  PFA,  PFD1 ,  PFD2 ,  R.  S,  P) 

ELSB  IF  (NATTYPB  .  BQ.  2  .AND.  NDFTYPB  . EQ.  1  .AND. 

1  ANSWER  .  BQ .  '  N '  )  THEN 

Sequential  attach  without  SLS  and  attach  at  target  unhnown 

CALL  SEQATKMAXR.  MAES,  PFA,  X.O,  PFD.  R,  S,  P) 

BLSB  IF  (NATTYPB  .  BQ.  2  .AND.  NDFTYPB  .BQ.  2  .AND. 

1  ANSWER  . BQ.  'N')  THEN 

Sequential  attach  with  SLS  and  attaoh  at  target  unhnown 

CALL  SBQATKMAER,  MAES,  PFA,  PFD1 ,  PFD2 ,  R,  S,  P) 

BLSB  IF  (NATTYPB  .BQ.  2  .AND.  NDFTYPE  .EQ.  2  .AND. 

1  ANSWER  .BQ.  'Y')  THEN 

Sequential  attaoh  with  SLS  and  attach  at  target  hnown 
CALL  8BQAT2(MAER,  MAES,  PFA,  PFD1 ,  PFD2 ,  R.  S.  P) 

BNDIP 

Specific  parameters  of  the  game 

WRITE(QPR, ' (A) ' )  'OThe  MINIMUM  attach  size  ?' 

RBAD( IN ,  FMT-* )  MINRV 

writb(QPR, '(A) ' )  'OThe  MAEIMUM  attach  size  ?' 

READ( IN ,  FMT-* )  MAERV 

WRITB(QPR, ' (A) ' )  'OThe  attach  Size  INCREMENT  ?' 


RSADCIN,  FMT- * )  INCRV 

WRITSCQPR, ' (A) ' )  'OThe  NUMBER  of  interceptors  ?' 

READ (IN,  fmt-*)  INT 

NRITSCQPR, ' (A) ' )  'OThe  TOTAL  NUMBER  of  targets  ?' 

RSADCIN,  FMT-*)  TARGETS 

WRITEC^PR, ' (A) ' )  'OThe  number  of  TYPES  of  targets  ?' 

RBADCIN,  FMT-*)  NTYPE 
c 

c  Only  one  target  type 

c 

IF  (NTYPB  . EQ.  I)  THEN 
VTYPBCl )-l 
NTAR(l) -TARGETS 
TVT-TARGETS 
c 

c  More  than  one  target  type 

o 

ELSE 

WRITBCQPR,  '(A)')  ' OBnter  first  the  RELATIVE  VALUE  and  the 
In  the  NUMBER  of  targets ' 

NRITBCQPR,  '(A)')  '  for  each  type.  Separate  the  two  entrie 
Is  for  eaoh  target  type  with' 

NRITB(§PR,  '(A)')  '  a  comma  and  hit  <CR>  following  the  ent 
lriee  for  eaoh  target  type: 
o 

o  Loop  through  each  target  type 

o 

DO  100  I  -  1,  NTYPE 

RBADCIN,  fmt-*)  VTYPB(I),  NT ARC I) 

TVT-TVT+ VTYPE ( I ) * HTAR C I ) 

TNT-TNT+NTARC I ) 

100  CONTINUE 

IF  (TNT  .NB.  TARGBTS)  THEN 

WRITBCTBR ,  '(A)’)  'OThe  sum  of  the  targets  in  the  indi 
lvldual  target  types  does  not ' 

NRITECTER,  '(A)')  '  equal  the  total  number  of  targets' 
STOP 
END  IF 
BNDIF 
o 

a  Calculate  the  fractional  values  for  VTYPB  and  NTAR 

o 

DO  110  I  -  1,  NTYPE 

VFRAC ( I ) - VTYPEC I ) * NTAR ( I ) / TVT 
NFRACC I ) -REAL  C  NTAR (I))/ TARGBTS 
110  CONTINUE 

a 

Q  Place  the  attach  sizes  in  an  array 

o 

A-l+C MAXRV-MINRV ) / INCRV 
DO  130  I-  1,  A 

RV ( I ) -MINRV+ ( I - 1 ) * INCRV 
120  CONTINUE 

0 


A- 6 


Call  Subroutine  BG  to  solve  for  the  basio  game 


CALL  BGCXBG.  YBG,  VBG,  B.  S,  P, 

1RV ,  A.  INT,  TARGETS . 

1  MAIN AT ,  MAXR,  MAIS ,  MAXNTYPE ,  NTYPE , 

1VFRAC.  NPRAC) 

Print  out  results  of  the  basic  game 

IP  (NOUT  . BQ.  1)  THEN 
O0T1-TBR 
OUT 2 -TER 

BLSE  IP  (NOUT  . EQ.  2)  THEN 
OCT 1 -PILE 
OOT2-FILE 

BLSE 

ODT1-TBR 

0UT2-FILE 

BNDIF 

Calculate  some  statistics  for  output 

DO  130  K-l.A 

DO  131  NT-1.  NTYPE 
TSU,NT)-0 

TMX-0 

TMY-0 

DO  132  1-0.  R 

DO  133  J-0,  S 

TS(K , NT)-XBG(K, NT, I)*P(I.J) *  YBG  (I,  NT , J)+ 

lTS(K.NT) 

IP  (I  .B9.  0)  THEN 

TMY-TMY+ J* YBG( K , NT , J ) 

BNDIF 

CONTINUE 

TMX-TMI+I *IBG(K .  NT , I) 

CONTINUE 

TS(t,NT)-TS(K,NT)*NTAR(NT) 

TMXBG ( I , NT ) -HINT C TMX • NTAR ( NT ) ) 

TMYBGC t . NT ) -HINT ( TMY ■ NTAR ( NT j ) 

CONTINUE 

CONTINUE 


DO  140  OUT -OUT 1 .  0UT2 

CALL  SUMMARY (NATTYPE , NDFTYPE , ANSWER , MINRV , KAXRV . INCRV , 
1  INT , TARGETS , R . S , PPA . PPD , PFD1 . PFD2 , OUT . 

1  MAINTYPB, NTYPE, VTYPE , NTAR) 

NAME- ‘ 1THE  ATTACIER ' ' S  BASIC  GAME  MINIMAI  STRATEGIES 
CALL  STPRINTCNAME , XBG , RV , A , R , OUT . MAIN AT , MAIR  . MAINTYPE , 
1  NTYPB.VFRAC, NPRAC) 


o  o  o 


NAME- ' 1THE  DEFENDER '' S  BASIC  GAME  MINIMAX  STRATEGIES 
CALL  STPRINT ( NAME , YBG , RV . A , S , OUT , MAXNAT , MAXS , MAXNTYPE , 

I  NTYPE , VFRAC , NFRAC ) 

o 

NAME-'ITHE  GAME  VALUES 
CALL  VPRINT ( NAME, VBG.RV, A, OUT, MAXNAT) 
c 

NAME-'ITHE  ATTACKER ' 'S  BASIC  GAME  TARGET  ALLOCATION  ' 

CALL  ALPRINTC NAME, XBG.RV, A, R, OUT, MAXNAT, MAIR, MAXNTYPE. 

1  NTYPE, VFRAC, NFRAC. TARGETS) 

c 

NAME- ' 1THB  DEFENDER  "  S  BASIC  GAME  TARGET  ALLOCATION  ' 

CALL  ALPRINT ( NAME . YBG , RV . A , S , OUT , MAXNAT , MAXS . MAXNTYPE , 

1  NTYPE , VFRAC . NFRAC , TARGETS ) 

c 

NAME-'ITHE  EXPECTED  NUMBER  OF  TARGETS  SURVIVING  ' 

CALL  ALVPRINT(NAME , VBG, RV, A, OUT, MAXNAT, 

1  MAXNTYPE. NTYPE, TS) 

o 

IF  (NTYPE  .NE.  1)  THEN 
ROB- . FALSE . 

NAME- ' 1THE  ATTACKER ' ' S  BASIC  GAME  RV  ALLOCATION  BY  TAR 

1GBT  TYPE' 

CALL  RVINTCOUNT ( NAME ,  TMXBG,  RV,  A.  OUT, 

1  MAXNAT , MAXNTYPE , NTYPE . ROB ) 

O 

WRITBCOUT, '(/)') 

NAME-'OTHE  DSFBNDSR ' ' S  BASIC  GAME  INTERCEPTOR  ALLOC AT I 
ION  BY  TARGET  TYPE' 

CALL  RVINTCOUNT (NAME.  TMYBG,  RV.  A,  OUT, 

I  MAXNAT,  MAXNTYPE.  NTYPE,  ROB) 

ENDIF 

140  CONTINUE 

Ask  for  the  number  of  robust  solutions  desired 

VRITE(QPR, ' (A) ' )  ' lPlease  enter  the  number  of  different  ra 
Ingee  of  RV'  s  for  which  robust  * 

WRITB(QPR, ' (A) ' )  '  solutions  are  to  be  found  ?’ 
WRITB(QPR, ' (A) ' )  '  *****  Enter  0  if  no  robust  solu 

Itlon  Is  desired** 

READ( IN , • )  J 
IF  (J  . BQ.  0)  THEN 
STOP 
ENDIF 

WRITB(QPR , ' (A) ' )  ' OThe  lower  and  upper  bounds  for  the  RV  range 
Is  aust  be  between' 

WRITB(QPR, ' (15,  "  and  '',15)’)  RV(1),  RV(A) 

DO  200  11-1.  J 

WRITE  (QPR,'(A)')  '0  The  lower  bound  :' 

READ  (IN,*)  MINRRV 

WRITE  (QPR , ' (A ) ' )  '0  The  upper  bound 
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READ  (IN,*)  MAXRRV 

AR- 1 + ( MAXRRV-MINRRV ) / INCRV 

DO  210  1-1.  AR 

210  RRV(I ) -MINRRV+ ( 1-1 ) * INCRV 

DO  220  1-1,  A 

IF  (MINRRV  . EQ.  RV(I))  THEN 
DO  221  12-1,  AR 

221  VBGRC 12 ) -VBG( I+I2-1 ) 

EKDIF 

220  CONTINUE 

c 

c  Find  the  robust  defender's  strategy,  YII 

c 

CALL  YROBUST(YII.  VBGR ,  R,  S.  P, 

1RRV,  AR.  I NT.  TARGETS. 

1MAXNAT . MAXR . MAIS , MAXNTYPE , NTYPE , NFRAC . VFRAC) 

c 

c  Find  the  Ill's  and  the  VII 's 

a 

CALL  XROBUST(XII,YII.  VII.R.S.P. 

1RV, A. INT, TARGETS. 

1MAXNAT . MAXR . MAIS . MAXNTYPE , NTYPE . NFRAC , VFRAC ) 

c 

DO  230  13-1. A 

DO  231  NT3-1 .  NTYPE 
TS(K3.NT3)-0 

TMX-0 

DO  232  13-0,  R 

DO  233  J3-0 .  S 

TS(K3 , NT3)-III(I3 , NT3 ,I3)*P(I3,J3)* 

1YII(NT3 . J3)+TS(K3 , NT3) 

233  CONTINUE 

TMX-TMX+I3*XII(K3 , NT3 ,13) 

232  CONTINUE 

TS ( 13 , NT3 ) -TS C 13 , NT3 ) * NTAR ( NT3 ) 

TMIII (13 . NT3 ) -NINT ( TMX * NTAR ( NT3 ) ) 

231  CONTINUE 

230  CONTINUE 

DO  240  NT3-1,  NTYPE 
TMY-0 

DO  241  J3-0.S 

TXY-TMY+J3*YII(NT3, J3) 

241  CONTINUE 

TMYII ( 1 . NT3 ) -NINT ( TMY • NTAR ( NT3 ) ) 


240  CONTINUE 

C 

o  Print  out  the  robust  solution 

c 

DO  250  OUT -OUT 1 .  OUT2 ,  FILE-TER 

0 

CALL  YRPRINTC YII ,  MINRRV,  MAIRRV ,  MAIS,  S,  OUT, 
1  MAZNTYPS,  NTYPE) 

C 


c 


o 


o 


a 


o 


o 

250 

200 


NAME- ' 1THE  OPTIMAL  ATTACK  STRATEGIES  AGAINST  THE  ROBUS 

IT  DEFENSE ' 

CALL  STPRINT ( NAME , III , RV . A , R , OUT , KAXNAT , MAXR , MAXNTYPE , 
1  NTYPE , VFRAC , NFRAC ) 

NAMB-'ITHB  EXPECTBD  TARGET  SURVIVAL  RATE  WITH  THE  ROBU 

1ST  DEFENSE' 

CALL  VPRINT(NAMB.VII,RV,A.0UT,MAXNAT) 

CALL  ALYRPRINT(YII,  MINRRV,  MAXRRV,  MAXS,  S,  OUT. 

1  MAXNTYPE,  NTT PB , NFRAC ,  TARGETS) 


NAME- ' 1THE  OPTIMAL  ATTACK  ALLOCATION  AGAINST  THE  ROBOS 

IT  DEFENSB ' 

CALL  ALPRINT  C  NAME , XII . RV , A , R , OUT . MAIN AT , MAXR , MAXNTYPE 
1  , NTYPE, VFRAC, NFRAC .TARGETS) 

NAME- ‘ 1THE  EXPECTED  NUMBER  OF  TARGETS  (SURVIVING)  WITH 
1  THE  ROBUST  DEFENSE' 

CALL  ALVPRINT ( NAME . VII , RV . A . OUT , KAXNAT . MAXNTYPE . NTYPE . 
1  TS) 

IF  (NTYPE  .NB.  1)  THEN 
ROB-. TRUE. 

NAME- ' 1THE  DBFENDER ' ' S  ROBUST  GAME  INTBRCEPTOR  ALL 
10CATI0N  BY  TARGBT  TYPE' 

CALL  RVINTCOUNT ( NAME ,  TMYII ,  RV.  A,  OUT. 

1  KAXNAT , MAXNTYPE , NTYPB . ROB ) 

WRITE (OUT,  '(/)') 


ROB- . FALSE . 

NAME-'OTHE  ATTACKER' '6  ROBUST  GAME  RV  ALLOCATION  B 

1Y  TARGBT  TYPE' 

CALL  RVINTCOUNT (NAME,  TMXII,  RV,  A.  OUT. 

1  KAXNAT, MAXNTYPE. NTYPE. ROB) 

END  IF 


CONTINUE 

CONTINUE 


STOP 

END 
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SUBROUTINE  SIMAT1(MAXR  , MAXS , PFA , PFD , R , S , P) 


Subroutine  SIMAT1  solves  for  the  probability  of  a  target 
surviving  under  simultaneous  attach  and  one  shot  defense 

INPUT  variables: 


MAXR 

INTEGER 

The  maximum  allowable  'R' 

MAXS 

INTEGER 

The  maximum  allowable  ’ S ' 

PFA 

REAL 

The  probability  that  a 

missile  will  fail  to  destroy 
a  target . 

PFD 

REAL 

The  probability  that  an 
interceptor  will  fail  to 
destroy  a  missile. 

R 

INTEGER 

The  maximum  number  of  missiles 
that  can  attach  a  target . 

S 

INTEGER 

The  maximum  number  of 

interceptors  that  can  attach 
a  target . 

OUTPUT  variable : 

P(I.J)  REAL  The  probability  that  a  target 

attached  by  I  missiles  and 
protected  by  J  Interceptors 
will  survive. 


INTEGER  R,  S 

REAL  P(0 : MAZE , 0 : MAXS) ,  MNFRA 

A-l-PFA 

D-l-PFD 

Loop  through  number  of  attaching  missiles 
DO  10  I  -  0,  R 

Loop  through  number  of  interceptors 

DO  20  J  -  0.  S 

Assign  value  to  P(I,J) 

IF  (I  . BQ.  0)  THEN 
P(I,J)  -  1 

ELSB  IF  (A  .HQ.  1.0  .AND.  D  .EQ.  1.0)  THEN 
IF  (I  .GT.  J)  THEN 
P(I.  J)  -  0 

ELSB 

P(I.  J)  -  1 

ENDIF 


ELSE 

MNINT  -  INTCJ/I) 


KNFRA  -  REAL(J)  'REALC I )-MNINT 
PI  -  (1-A*(1.0-D)**(MNINT+1))**(I*MNFRA) 
P2  -  (1-A*(1.0-D)*  *MNINT) *  * ( I* ( 1-MNPRA) ) 
P(I.J)  -  PI  *  P2 
SNDIF 

CONTINUB 

CONTINUE 


RETURN 


SUBROUTINE  SIMAT2( MAXR , MAXS , PFA , PFD1 , PFD2 , 
R.S.P) 


Subroutine  SIMAT2  calculates  for  the  probability  of  a  target 
surviving  under  simultaneous  attack  and  shoot,  look,  shoot 
defense . 

INPUT  variables : 


MAXR 

INTEGER 

The  maximum  allowable  'R' 

MAXS 

INTEGER 

The  maximum  allowable  ' S ’ 

PFA 

REAL 

The  probability  that  a 
missile  will  fail  to 
destroy  a  target. 

PFD1 

REAL 

The  probability  that  a  first 
salvo  interceptor  will  fail 
to  destroy  a  missile 

PFD2 

REAL 

The  probability  that  a  second 
salvo  interceptor  will  fail 
to  destroy  a  missile 

R 

INTEGER 

The  maximum  number  of 

missiles  that  can  attack  a 
target 

S 

INTEGER 

The  maximum  number  of 
interceptors  that  can 
defend  a  target 

OUTPUT  variables : 

P(I,J) 

REAL 

The  probability  that  a  target 

attacked  by  I  missiles  and 
protected  by  J  Interceptors 
will  survive 


INTEGER  R , S 

REAL  P( 0 : MAXR , 0 : MAXS ) 

DEFINITION  OF  LOCAL  VARIABLES: 

INPUTS  ARE  (1)  MAX  NUMBER  OF 

ATTACKERS  N.  (2)  MAX  NO.  OF  INTERCEPTORS  M,  (3)  MAX  NO. 

OF  SALVOS  SAL,  (4)  SINGLE-SHOT  INTERCEPTOR  KILL  PROBABILITY 

ON  EACH  SALVO  PKILL(l) . PKILL(SAL),  AND  (5)  ATTACKER 

RELIABILITY  RHO.  OUTPUTS  ARE  (1)  MINIMUM  PROBABILITY  OF 
TARGET  DESTRUCTION  WHEN  K  SALVOS  REMAIN,  S2(I,J), 

1-1 . M;  J-l . N;  K-I . SAL,  (2)  OPTIMAL  NUMBER 

OF  INTERCEPTORS  TO  USE  UNDER  EACH  CONDITION  D(I.J),  AND 
(3)  EXPECTED  NUMBER  OF  INTERCEPTORS  REMAINING  AT  THE 
END  OF  THE  ATTACK  W2(I,J).  IN  CASE  OF  TIES,  THE 
CHOICE  IS  MADE  TO  MAXIMIZE  V2(I,J). 

ARRAYS  NEEDED  ARE  S1(M+1,N+1) ,S2(M+1 ,  N+ 1 ) ,  D(M-1,N*1). 

PA(M*1 , N+l ) ,  RRCN+l , N) ,  TOCN),  TI(N).  PKILL(R) , 


C  W1 (M* 1 , N+l ) ,  AND  W2(M+1 . N+l ) . 

C 

DIMENSION  S1(51,51),S2(51.51).D(51.51).W1(51,51).W2(51.51). 
1  PA(51,51).RR(51,50) , T0( 30) , T1 ( 50) , PKILLC 5 ) 

INTEGER  D, SAL 
C 

C  SET  LOCAL  INPOT  VARIABLES  TO  THEIR  CORRESPONDING  'SEARED'  INPOT 

C  VARIABLES 

C 

N-R 
M-S 
SAL-2 

RHO-l.-PFA 
PKILLC 1 )-l , -PPD1 
PKILL(2)-1 . -PFD2 
C 

C  PKILLC 1)  IS  THE  S-S  KILL 

C  PKILLC 2)  IS  THE  S-S  KILL 

C 

C  INITIALIZE;  DO  CASE  OF  1 
C 

RHOBAR  -  1.  -  RHO 
Ml-M+1 
Nl-N-i-1 

DO  200  I-1.M1 
W2(I,1)  -  1-1 
200  S2CI.1)  -  0. 

S2( 1 , 2) -RHO 
DO  210  J-3.N1 
Jl-J-1 

210  S2(l, J)-l. -RHOBAR* (1.-S2C1. 01)) 

DO  211  I-1.M1 
DO  211  J-2.N1 
W2CI.J)  -  0.0 

211  CONTINUE 
9-1 . -PKILLC SAL) 

DO  220  J-2.N1 
IF( J  .GE.  2)  GO  TO  1211 
DO  1200  1-2.  Ml 
D(I.J)  -  1-1 

1200  S2( 1 , 2 )  -  9**(I-1) 

GO  TO  220 
1211  Jl-J-1 
ZJ1-J1 

DO  220  1-2, Ml 
Il-I-l 
ZI1-I1 
ZIJ-ZI1/ZJ1 
INTIJ-I1/J1 
Z INTI J- INTI J*1 . E-2 
IP  CZINTIJ  .LT.  ZIJ)  GO  TO  212 
C 

C  ZIJ  IS  AN  INTEGER 
C 


PROB.  ON  THE  FIRST  SALVO . 

PROB.  ON  THE  SECOND  SALVO. 

SALVO  LEFT. 


S2(I.J)-1.-(1.-RH0*(Q,*INTIJ;;**J1 
GO  TO  214 
212  CONTINUE 

ZIJ  IS  NOT  AN  INTEGER 

INTI J1 -INTI J* 1 
IEXP1-I1-J1 *  INTI J 
IEXP2-J1 -IEXP1 

FI -Cl .  -RH0*CQ“INTIJ1))**IEXP1 
F2-( 1 . -RSO* (Q*  * INTIJ ) ) *  * IEXP2 
S2( I , J)-l • -FI  *  F2 
214  CONTINUE 
D(I.J)-H 
220  CONTINUE 
K-l 

NOW  WRITE  RESULTS  FOR  X-l 
230  CONTINUE 

IF  (SAL  . £Q .  1)  GO  TO  900 

MAIN  LOOP  FOR  K  SALVOES.  K>1 

DO  420  K-2.SAL 
XX  -  SAL-Kfl 
g-l.-PKILL(KK) 

DO  300  I-1.M1 
DO  300  J-1.N1 
Wl(I.J)  -  W2(I.J) 

300  S1(I. J)  -  S2CI.J) 

FIRST  FIX  J 

DO  400  J-2.N1 

Jl-J-1 

ZJ1-J1 

CALCULATE  ALL  PA(-.-)  VALUES 

PA( 1 . J)-l • 

DO  310  JJ-1 , J1 
310  PA( 1 , JJ ) -0 . 

PIRST  CALCULATE  TO  AND  T1  VECTORS 

DO  380  1-2. Ml 

Il-I-l 

ZI1-I1 

ZIJ-ZI1/ZJ1 

INTIJ-I1/J1 

ZINTIJ-INTIJ-*- 1 .  E-5 

IF  (ZINTIJ  .LT.  ZIJ)  GO  TO  330 


ZIJ  IS  AN  INTEGER 

FAC-g* ’INTIJ 
DO  320  L-1.J1 
T1(L)-FAC 
320  T0(L)-1 . -FAC 
GO  TO  335 
330  CONTINUE 

ZIJ  IS  NOT  AN  INTEGER 

LIM-J1-I1+ J1 • INTIJ 
LIM1-LIM+1 
FAC-g** INTIJ 
DO  332  L-l.LIM 
T1(L)-FAC 

332  T0(L)-1.-FAC 
FAC-Q*FAC 

DO  333  L-LIM1.J1 
T1(L)-FAC 

333  T0(L)-1 . -FAC 
335  CONTINUE 

NOW  FINISHED  WITH  TO  AND  T1 

RR(l.l)-TOU) 

RR(2.1)-T1(1) 

IF  (J1  -Eg.  X)  GO  TO  360 

DO  340  L-2.J1 

Ll-L+1 

LM1-L-1 

RR(1.L)-RR(1.LM1)*T0(L) 

340  RR(L1 ,L)-RRCL,LM1)*T1(L) 

DO  350  L-2.J1 
LM1-L-1 
DO  350  J2-2.L 
JM1-J2-X 

RR(J2.L)«RR(J2.LM1)*T0(L)+RR(JM1 ,LM1)*T1(L) 

350  CONTINUE 
360  CONTINUE 

DO  370  JJ-l.J 
370  PA(I. JJ)-RR(JJ. Jl) 

360  CONTINUE 

ALL  PA( - , - )  VALUES  ARE  NOW  READ* 

TEST  IS  TENTATIVE  VALUE  OF  S(I,J) 

ITEST  IS  TENTATIVE  VALUE  OF  D(I,J) 

WTEST  IS  TENTATIVE  VALUE  OF  W(I,J) 

382  CONTINUE 

DO  395  1-2, Ml 
TEST-1 . 


rS-A 


DO  390  II-l.I 
MI-I-II+1 
PAG  -  0.0 
DO  385  JJ-l.J 

385  FAC-FAC+PACII. JJ)*S1(MI. JJ) 
TEST-AMINlC FAC. TEST) 

£ 

C  UPDATE  TENTATIVE  VALUES 

C 

390  CONTINUE 

S2( I , J) -TEST 
395  CONTINUE 
400  CONTINUE 
C 

c  NOW  WRITE  RESULTS 

C 

420  CONTINUE 
C 

500  CONTINUE 
C 

c  STORE  PIJ'S 

C 

DO  154  I-O.N 
DO  155  J-O.M 

P(I. J)-1-S2(J+1.I-1) 
155  CONTINUE 
154  CONTINUE 

RETURN 
900  END 


yy 


SUBROUTINE  SB9AT1(MAXR , MAXS , PFA . PFD1 , PFD2 , R , S , P) 

Subroutine  SBQAT1  solves  for  the  probability  of  a  target 
surviving  under  unknovn  sequential  attaok  and  one  shot  or  shoot , 
look,  shoot  defense 

INPUT  vadriables: 


o 

MAZR 

INTEGER 

The  maximum  allowable  'R' 

o 

MAZS 

INTBGBR 

The  maximum  allowable  ' S ' 

o 

PFA 

REAL 

The  probability  that  a 

0 

missile  will  fall  to 

o 

destroy  a  target 

0 

PFD1 

REAL 

The  probability  that  a 

o 

first  salvo  interceptor 

0 

will  fall  to  destroy  a 

o 

missile 

o 

PFD2 

REAL 

The  probability  that  a 

0 

seoond  salvo  lnteroeptor 

o 

will  fail  to  destroy  a 

o 

missile 

0 

R 

INTBGBR 

The  maximum  number  of 

o 

missiles  that  can  attack 

0 

a  target 

0 

S 

INTEGER 

The  maximum  number  of 

0 

Interoeptors  that  may 

o 

Q 

attaok  a  target 

o 

OUTPUT  variable: 

o 

o 

P(I.J) 

REAL 

The  probability  that  a 

target  attaoked  by  I 
Bissiles  and  proteoted 
by  J  interoeptors  will 
survive 


INTBGBR  R.S 

REAL  PCO.MAXR,  O.MAXS) 

LOCAL  VARIABLE  DEFINITIONS 

RR( . . D)  IS  THB  BXPBCTBD  DAMAGE  CURVE  WITH  ONE  SHOT  LEFT 
9(. ,D)  IS  THB  BZPBCTBD  DAMAGE  CURVE  WITH  TWO  SHOTS  LEFT 
SLOPE (D)  YIELDS  PRIM  READ  SLOPS 

NQ(D)  IS  THB  NUMBER  OF  INT  TO  SHOOT  AT  1ST  RV.  1ST  VOLLEY 
NR(D)  IS  THB  NUMBER  OF  INT  TO  SHOOT  AT  NEXT  RV,  2ND  VOLLEY 
NRR(D)  IS  THB  NUMBER  REMAINING  IF  MISS  (D-NQ(D) ) 

DIMENSION  9(0 : 100 , 0 : 100 ) , RRCO : 100 . 0 : 100 ) , SLOPS( 0 : 100) , 

1  NR(0 : 100)  ,  NQ(0  •'  100)  . RTRY(0  : 100 )  ,  9TRY(0 : 100)  , 

1  NRR(0:100) 


CAPBSS-1 . -PFA 
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o  o  o 


9(IA,ID)-RR(IA.ID) 

SLPTRX-Q(IA.ID)/IA 
IF (SLPTRY.GT. SLPMAX)  SLPMAX-SLPTRY 
END  DO 
NQ(ID)-0 
DO  IK-1,  ID 
SLPMAX-0 . 

POW-PFD1 *  *  IK 
DO  IA-l.R 

9TRY(IA)-POW*RR(IA, ID-IK)+( 1 . -PON) *Q(IA-1 . ID- IK) 
SLPTRX-QTRY ( I A ) / IA 
IF ( SLPTRY . GT . SLPMAX )  SLPMAX-SLPTRY 
END  DO 

IF ( SLPMAX. LT.SLOPBC ID))  THEN 
DO  L-l.R 

Q(L. ID)-QTRY(L) 

END  DO 

SLOPE (ID) -SLPMAX 
NQ(ID)-IK 
END  IF 
END  DO 
END  DO 

RETURN  THE  SURVIVAL  PROBABILITIES 

DO  I-O.R 

DO  J-O.S 

P(I,J)-1.-Q(I,J) 

END  DO 
END  DO 

DO  ID-1 . S 

NRRCID)-ID-NQ(ID) 

END  DO 

RETURN 

END 
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o  o 


SUBROUTINE  SEQAT2 ( MAXR , MAXS ,  PF A .  PFD 1 ,  PFD2  ,  R ,  S  ,  P  ) 

C 

o  Subroutine  SEQAT2  solves  for  the  probability  of  a  target 

o  surviving  under  knovn  sequential  attack  and  shoot,  look,  shoot, 

o  defense 

c 

c  INPUT  variables: 

c 


c 

MAXR 

INTEGER 

The  maximum  allowable  'R' 

0 

MAXS 

INTEGER 

The  maximum  allowable  ' S ' 

c 

PFA 

REAL 

The  probability  that  a 

c 

missile  will  fail  to 

c 

destroy  a  target 

c 

PFD1 

REAL 

The  probability  that  a 

c 

first  salvo  interceptor 

a 

will  fail  to  destroy  a 

c 

missile 

c 

PFD2 

REAL 

The  probability  that  a 

0 

second  salvo  Interceptor 

c 

will  fall  to  destroy  a 

o 

missile 

c 

R 

INTEGER 

The  maximum  number  of 

0 

missiles  that  can  attack 

0 

a  target 

0 

S 

INTEGER 

The  maximum  number  of 

0 

interceptors  that  can 

a 

A 

defend  a  target 

W 

o 

OUTPUT  variable: 

o 

0 

P(I. J) 

REAL 

The  probability  that  a 

0 

target  attacked  by  I 

o 

mlSBlles  and  protected 

o 

by  J  interceptors  will 

0 

o 

survive 

c 

INTEGER  R.S 

REAL  P  (  0  •'  MAXR ,  0 

:  MAXS ) 

c 

O  LOCAL  PROGRAM  NOTES: 

O 

O  COMPUTES  AN  OPTIMAL  PRIM  READ  DEFENSE  WITH  A 

O  SLS  CAPABILITY  UNDER  THE  ASSUMPTION  THAT  THE 

O  DBFBNSB  CAN  DETERMINE  THE  SIZE  OF  THE  ATTACK 

O  AT  THE  MOMENT  THAT  THE  ATTACK  BEGINS. 

O 

O 

O  DBS( .  ,  .  )  IS  THE  EXPECTED  DESTRUCTION 

o 

REAL  SUR( 100 . 100 ) . DES( 100 . 100) . J RSLPBC 100 ) 
INTBGER  DBEC100, 100) ,EEE(100. 100) 

PARAMETER  EPS I  IS  A  TOLERANCE 


o 


EPS I -0.000001 


SLOPE-1. /R 

IIPRIM  WILL  REMAIN  1  WHEN  DBS < -JA* SLOPE  FOR  ALL  JA 

IIPRIM-1 
DO  ID-l.S+1 
S0R(1,ID)-1. 

DBS( 1 , ID)-0 . 

UO  DO  JA-2.R+1 

SUR( JA, ID)-PFA*SURC JA-1 . ID) 

DESC JA, ID)-1 . -SURC JA, ID) 

IPRIM  WILL  BE  SET  TO  1  IP  PRIMRBAD  OK  AT  JA.ID 


IPRIM-0 
DEBC JA, ID)-0 
EBEC JA, ID)-0 
SURTRY-SURC JA.ID) 

DO  LODES- 1. ID 

DO  LOBBE-l.ID-LODEE+1 

INDTEM-ID-LODBE-LOEBE+2 

TRYSUR- ( 1 . -PFD1 • • ( LODES- 1 ) ) • SUR( JA- 1 , ID-LODEE+ 1 ) + 
C  PFD1  *  • CLODEE-1 )  )  • 

( 1 . -CPFD2*  *  CLOBBE-1) ) * ( 1 . -PPA) ) * 

SURC JA-1, INDTBM) 


O  TBST  FOR  UPDATE:  “THRIFTY  ALLOCATION"  BREAKS  TIES 
C  IN  FAVOR  OF  THE  SMALLEST  POTENTIAL 

O  ALLOCATIONS.  IF  STILL  TIED,  IT  PICKS 

O  THE  SMALLEST  INITIAL  VOLLEY, 

o 

SRTRYN-SURTRY-BPSI 

IF ( TRYSUR. LT.SRTRYN)  GO  TO  1000 

SRTRYP-SURTRY+EPSI 

IF ( TRYSUR. LE.SRTRYP)  THEN 

ISMINC-DEE(JA,ID)+EBE(JA,ID)+2 
ISMTRY-LODEB+LOBBE 
IFCISMINC.LT. I SMTRY)  GO  TO  1000 
IF ( ISMINC . BQ . I SMTRY )  THEN 

IFCDEBC JA, ID) . LE . LODBE-1)  GO  TO  1000 
BND  IF 
END  IF 

C  UPDATE  IS  NBCBS8ARY 

SUR( JA.ID) -TRYSUR 
DESC JA , ID)-1 . -SURC JA, ID) 

DBEC JA. ID) -LODES- 1 
BBBCJA, ID)-LOEB8-l 
SURTRY-TRYSUR 

C  UPDATB  UNECBSSARY 
1000  CONTINUE 

BND  DO 


END  DO 
BND  DO 
END  DO 
o 

O  DETERMINE  WORST  SLOPE  FOR  EACH  DEFENSE  SIZE 

o 

DO  ID-1, S+l 

PRSLPE(ID)-0. 

DO  IA-2.R+1 

SLPTRY-DESC IA , ID)/(IA 

IF  C  SLPTRY . GT . PRSLPE ( ID ) )  PRSLPE ( ID ) - SLPTRY 
END  DO 
END  DO 
c 

O  RETURN  PI J ' S 
O 

DO  I-O.R 
DO  J-O.S 

P(I, J)-1-DES(I+1, J+l) 

END  DO 
BND  DO 
O 

RETURN 

BND 


APPENDIX  C 


SUBROUTINES  BG,  YROBUST,  AND  XROBUST 


fhX 


c 


SUBROUTINE  BGCZBG,  YBG ,  VBG ,  R,  S,  P, 

1RV,  A.  INT,  TARGETS, 

1MAT,  MAZR,  MASS,  MAZNTYPB .  NTYPE , 

1VFRAC,  NFRAC) 

C 

c  Subroutine  BG  finds  and  returns  the  solution(s)  for  the  basic  game 

c 

c. . . . . . . 


c 

c  INPUT  variables : 


c 


c 

R 

INTEGER 

Maximum  RV's  at  a  single 

G 

target 

c 

s 

INTEGER 

Maximum  interceptors  at  a 

c 

single  target 

c 

P(I.J) 

REAL 

Probability  of  survival  for 

c 

a  target  attacked  by  I 

RV ' s  and  protected  by  J 

o 

0 

interceptors 

o 

RV 

INTEGER 

The  attack  sizes 

c 

A 

INTEGER 

The  number  of  attack  sizes 

c 

INT 

INTEGER 

Number  of  interceptors 

c 

TARGETS 

INTEGER 

Number  of  targets 

c 

MAT 

INTEGER 

Maximum  permissable  A 

c 

MAZR 

INTEGER 

Maximum  permissable  R 

c 

MAZS 

INTEGER 

Maximum  permissable  S 

c 

MAZNTYPB 

INTEGER 

Maximum  permissable  NTYPE 

c 

NTYPE 

INTEGER 

The  number  of  target  types 

o 

VFRAC(Z) 

REAL 

Fractional  value  of  all 

o 

NFRAC ( K ) 

type  K  targets 

o 

REAL 

Fractional  number  of  all 

o 

o 

- 

type  Z  targets 

0 

0 

Input  variable 

type  declaration 

c 

INTEGER  R.  S,  A,  INT,  TARGETS 
INTEGER  MAT.  MAXR,  MAZS 
INTEGER  RV(MAT) ,  MAZNTYPB ,  NTYPE 
REAL  P  C  0 : KAZR , 0 : MAZS ) 

REAL  VFRAC(MAZNTYPB) .  NFRAC (MAZNTYPB) 

C 

. . . . . 

o 

o  OUTPUT  variables : 

c 

o  ZBG(A,  K,  I) 

c 
o 
c 

O  YBG(A,  Z,  J) 

C 


REAL  Min-max  strategy  for  the 

attacker  for  attack  size 
A:  fraction  of  K'th  type 
targets  with  I  RV's 

REAL  Min-max  strategy  for  the 

defender  for  attack  size 


A:  fraction  of  K'th  type 
targets  with  J  interceptors 

VBG(A)  REAL  Ga me  value  associated  with 

X3G  and  YBG 


Output  variable  type  declaration 

REAL  XBG(MAT ,  MAXNTYPE ,  0  : MAXR )  ,  YBG (MAT, MAXNTYPE, O-MAXS) 
REAL  VBG(MAT) 


Local  variables  needed  for  XMP .  See  XMP  Dictionary  for 
definitions . 


Integer  parameters 

INTEGER  MAXM , MAXN , MAXA , COLMAX . PICE 

PARAMETER  (MAXM-330,  MAXN-1000,  MAZA-30000.  COLMAJ-550, 
1PICK-7) 

INTEGER  LENI1,  LENMI1,  LENMR1 ,  LENR1 

PARAMETER  ( LENI 1 -60000 . LENMI 1 -9 . LENMR1-8 , LENR1 -60000 ) 

Double  Precision  arrays  and  variables 

DOUBLE  PRECISION  B(MAXM) , BASCB(MAXM) , BASLB(MAXM) , BASUB(MAXM) 
DOUBLE  PRECISION  BLOW(l ), BOUND, CANDA (COLMAX . PICK) , CANDCJC PICK) 
DOUBLE  PRECISION  CJ . COLA(COLMAX) , LJ , MEKR(LENR1 ) , UJ , UZERO(MAXM) 
DOUBLE  PRECISION  XBZ ERO( MAXM ) ,YQ (MAXM) , Z 

Integer  arrays  and  variables 

INTEGER  BNDTYP , COLLEN , DFEASQ , DTERM , DUMBR , ERROR , FACTOR 
INTEGBR  IOERR , IOLOG , ITER . ITER 1 , ITERS , LOOK , M , MAPI ( LENMI 1 ) 
INTEGER  MAPR(LENMRl) .N.NTYPES, PRINT, TERMIN 
INTEGER  UNBDDQ,MEMI(LENI1) , BASIS (MAXM) , CAND(PICK) 

INTEGER  CANDI ( COLMAX , PICK ) , CANDL( PICK ) , COLI ( COLMAX ) 

INTBGBR  ROWTYP(MAXM) ,  STATUS (MAIN) 

I/O  numbers  for  input,  error  output,  and  results  output 

COMMON/ IO.UNIT/ IOIN , IOERR , IOLOG 

Other  ooamon  variables 


COMMON  /XMPLEN/  LENI.  LENMI,  LENMR,  LENR 
INTEGER  LENI,  LENMI,  LENMR,  LENR 


LENI-L3NI1 

LENMI-LENMI1 

LENR-LENR1 

LENMR-LENHR1 

LOOK-200 

FACTOR-60 

PRINT-0 

BNDTYP-1 

CALL  XMAPS ( BNDTYP . IOERR , MAPI , MAPR , MAXA , MAXM , MAYN , MEMI , MEMR ) 

Prepare  IMP  to  solve  for  the  first  attack  size  RV(l) 

Set  row  type  and  right  hand  side  for  the  inequality  constraints 

Nl-NTYPE* (R+l ) 

DO  100  1-1,  NX 
B(I)-0 

RONTYP(I)-+l 

CONTINUE 

Set  row  type  and  right  hand  side  for  the  equality  constraints 

N2-NTYPB* (R+2) 

DO  110  I-Nl+1 ,  N2 
B(I)-1 
ROWTYP(I)-0 
CONTINUE 

B ( N2+ 1 ) -REAL ( I NT ) / REAL ( TARGETS ) 

ROWTYP(N2+1)-0 

Let  N  be  the  current  number  of  variables;  Incremented  by  XADDAJ 
N-0 

Construct  the  first  NTYFE*(S+1)  structural  variables  associated 
with  the  defender's  strategy 

DO  120  NT-1.  NTYPE 
DO  121  J-0.  S 
COLLEN-O 
DO  122  1-0,  R 

IP  (P(I.J)  .NE.  0)  THEN 
COLLEN -COLLBN + 1 

COLA ( COLLEN )  —  P( I . J ) * VFRAC  (  NT  ) 

COLI ( COLLEN ) -1+ 1 +( R+ 1 )*( NT- 1 ) 

ENDIF 
CONTINUE 
COLLEN - COLLEN + 1 


Set  the  number  of  constraints 


M-N2+1 

Construct  the  slack,  surplus  and  artificial  variables  with  XSLACK 

CALL  XSLACK (B , BASCB , BASIS , BASLB , BASUB , BLOW . 

BNDTYP ,  BOUND , 

COLA , COL I . COLMAX , IOERR , 

M ,  MAPI , MAPR , MAXM , MAXN , MEMI , MEMR . 

N ,  ROWTYP , STATUS , UZERO , XBZERO , Z ) 

Solve  the  LP  with  the  primal  simplex  method 

CALL  XPRIML(B , BASCB , BASIS , BASLB , BASUB . BNDTYP , BOUND , 

CAND , CANDA , CANDC J , CANDI , CANDL , 

COLA . COL I , COLMAX , 

PACTOR , IOERR , BOLOG . ITER1 , ITER2 , 

LOOK , M , MAPI , MAPR , MAXM , MAZN , MEMI . MEMR , 

N , NTXPE2 . PICK , PRINT . STATUS , TERMIN , UNBDDQ , 

0ZERO, XBZERO, VQ.Z) 

Check  if  solution  has  been  found 

IF  (TERMIN  .NE.  1)  THEN 

WRITE  (IOERR, '(A)')  ' 1 THERE  IS  AN  ERROR  IN  THE  FORMULATION 

WRITE  (IOERR, '(A. 14)')  '  TERMINATION  CODE  -  ',  TERMIN 
STOP 
ENDIF 

Store  in  ZBG,  YBG,  and  VBG  the  solutions  to  the  LP  for  RV(i) 

DO  200  NT-1,  NTTPE 
DO  201  I-  0,  R 

XBG( 1 , NT, I)-UZERO(I+l+(NT-l ) * (R+l ) ) 

CONTINUE 

CONTINUE 

DO  210  I-  1.  N2+1 

IF  (BASIS(I)  .LE.  ( S+l ) *NTYPE)  THEN 
YBG(1,(BASIS(I)-1)/(S+1)+1. 

M0D(BASIS(I)-1, (S+l)))-XBZERO(I) 

ENDIF 

CONTINUE 

VBG( 1 )-Z 

Solve  for  the  remaining  RV's 

IF  (A  .Eq.  1)  THEN 
RETURN 
ENDIF 

DO  300  I  -  2,  A 


C-4 


COLA (COLLEN)- 1.0 
COLIC COLLEN )-Nl+ NT 
IP  (J  . NB .  0)  THEN 
COLLEN -COLLEN+1 
COLA  ( COLLEN ) -J  *  NFRAC ( NT ) 

COLI ( COLLEN )-N2+l 
ENDIF 

Call  XADDAJ  to  enter  the  column  for  YJ  Cor  variable  J+l) 
OBJ-O 

CALL  XADDAJ(OBJ . COLA , COLI , COLLEN . COLMAX , IOERR , 

X , MAPI . MAPR , MEMI , MEMR , N) 

CONTINUE 

CONTINUE 

Construct  the  columns  for  the  ' s '  variables 

DO  130  NT-1,  NTYPE 
DO  131  I-l.R+1 
COLACD-l 

COLI(I)-I+(NT-l)*(R+l) 

CONTINUE 

Call  XADDAJ  to  enter  the  column  for  s  Cor  variable  S+2) 

COLLEN-R+1 

OBJ-1 

CALL  XADDAJ  C  OBJ . COLA , COLI , COLLEN , COLMAX . IOERR , 

X . MAPI . MAPR . MEMI . MEMR . N ) 

CONTINUE 

Construct  the  column  for  the  't'  variable 
COLLEN-O 

DO  140  NT-1,  NTYPE 
DO  141  I-l.R 

COLLEN -COLLEN  + 1 

COLA  C  COLLEN  )— I  •  NFRAC  C  NT ) 

COLI ( COLLEN )-I+l+( NT- 1 ) * (R+l ) 

CONTINUE 

CONTINUE 

Call  XADDAJ  to  enter  the  column  for  t  (or  variable  S+3) 

OBJ-RBALC -RV( 1 ) ) /REAL( TARGETS ) 

CALL  XADDAJ ( OBJ , COLA , COLI , COLLEN , COLMAX , IOERR , 

X . MAPI . MAPR . MEMI . MEMR , N) 

Start  all  the  structural  variables  at  their  lower  bounds 

DO  150  J-l,  N 
STATUS( J)-0 


J-NTYPE* (S+2)+l 

CJ-REALC -RV( I ) ) /REAlf TARGETS ) 

CALL  ZCEGCJCCJ , J . MAPI , KAPR , MEMI , MEMR ) 


C 


310 

O 


c 

c 

o 


321 

320 


330 


DO  310  11-1, M 
J-BASIS(Il) 

CALL  XGETA J ( C J , COLA , COLI , COLLEN , COLMAX , 
1  IOERR, J, MAPI. MAPR, MEMI, MEMR) 

BASCB(I1)-CJ 

CONTINUE 


1 

1 

1 

1 

1 

1 


CALL  XPRIMLCB , BASCB . BASIS . BASLB , BASUB . BNDTYP , BOUND , 
CAND , CANDA , CANDCJ , CANDI . CANDL , 

COLA, COLI, COLMAX, 

PACTOR . IOERR , EOLOG . ITER1 , ITER2 . 

LOOK . M , MAPI , MAPR . MAXM . MAXN . MEMI . MEMR . 

N , NTYPE2 , PICX , PRINT , STATUS , TBRMIN , UNBDDQ , 
UZERO , XBZERO , YQ , Z ) 


Store  In  XBG.  YBG,  and  VBG  the  solutions  to  the  LP  for  RV(i) 


DO  320  NT-1,  NTYPB 
DO  321  12-  0,  R 

XBG(I,NT,I2)-UZBR0(I2+D(HT-1)*(R+1)) 

CONTINUE 

CONTINUE 

DO  330  12-  1,  N2+1 

IP  (BASISU2)  .LE.  (S+1)*NTYPE)  THEN 
YBG(I ,  (BASIS(I2)-1)/(S+1)+1 , 

1M0D( BASIS(I2)-1 , ( S+ 1 ) ) )-XBZERO( 12) 

ENDIF 

CONTINUE 


VBGCD-Z 

300  CONTINUE 

O 

o 

RETURN 

O 

END 


O 

o  Subroutines  for  changing  an  objective  row  coefficient 

o 


SUBROUTINE  XCHGCJC CJ , J , MAPI , MAPR , MEMI , MEMR ) 
COMMON / XMPLBN / LEN I , LENMI , LENMR , LENR 
DOUBLE  PRECISION  C J , MEMR ( LENR ) 

INTEGER  MAPI C LENMI ) ,  MAPR(LENMR) 

INTEGER  MEMI  (LEND 


C-6 


SUBROUTINE  YROBUST(YII,  VBG ,  R.  S,  P, 

1RV ,  A,  1ST,  TARGETS, 

1MAT,  KAXR ,  MAYS,  MAYNTYPE,  NTYPE ,  NFRAC,  VPRAC) 

C 

o  Subroutine  ROBUST  finds  and  returns  the  robust  defender's  strategy 

c  over  a  range  of  attach  sizes 

c 

C*************************  PRIMAL  VERSION  ****#******«##*####«**#*#*#* 


INPUT  variables : 

VBG(A) 

REAL 

Game  value  associated  with 

YBG  and  YBG 

R 

INTEGER 

Maximum  RV ' s  at  a  single 
target 

S 

INTBGER 

Maximum  interceptors  at  a 
single  target 

P(I. J) 

REAL 

Probability  of  survival  for 
a  target  attacked  by  I  RV ' s 
and  protected  by  3 
interoeptors 

RV 

INTEGER 

The  attack  sizes  over  which 
the  robust  defense  will  be 
defined 

A 

INTEGER 

The  number  of  attack  sizes 

INT 

INTBGER 

Number  of  interceptors 

TARGETS 

INTBGER 

Number  of  targets 

MAT 

INTEGER 

Maximum  permissable  A 

MAYR 

INTEGER 

Maximum  permissable  R 

MAYS 

INTEGER 

Maximum  permissable  S 

MAYNTYPE 

INTEGER 

Maximum  permissable  NTYPE 

NTYPE 

INTEGER 

Number  of  target  types 

NFRAC(K) 

REAL 

The  fraction  of  targets  of 
type  Y 

VFRAC(Y) 

REAL 

The  fraction  of  the  total 
value  for  type  Y  targets 

Input  variable  type  declaration 

INTEGER  R,  S,  A,  1ST,  TARGETS 

INTEGER  MAT,  MAYR,  MAYS,  MAYNTYPE,  NTYPE 

INTEGER  RV(MAT) 

REAL  P(0 : MAYR , 0 : MAYS ) ,  VBG (MAT) ,  NFRACC MAYNTYPE) 
REAL  VPRACC MAYNTYPE) 


OUTPUT  variables 


o 


ooo  o  o  o 


YII(K.J) 


REAL 


I 


c 

v 


•• 

1 


I 


o 

0 

o 

0 

O 

O 

O 

c 

o 

o 

o 

o 

O' 

c 

c 

o 

o 

o 


Robust  strategy  for  the 
defender  over  attach  sizes 
RV:  fraction  of  type  K 
targets  with  J 
interceptors 


Output  variable  type  declaration 
REAL  YlKMAXNTYPE  ,  0  :  MAXS) 


Local  variables  needed  for  XMP .  See  XMP  Dictionary  for 
definitions. 


Integer  parameters 

INTEGER  MAXM , MAXN , MAXA , COLMAX , PICE 
PARAMETER  (MAXM-1500,  MAXN-1700,  MAXA-20000) 

PARAMETER  (COLMAX-lOO,  PICK-7) 

INTEGER  LENI1,  LENMI1,  LBNMR1 .  LENR1 

PARAMETER  ( LENI 1-100000, LBNMI 1 -9 , LENMR 1 -8 , LENR 1 - 80000 ) 


Double  Precision  arrays  and  variables 

DOUBLE  PRECISION  B(MAXM) , BASCB(MAIM) , BASLB(MAXM) , BASUB(MAXM) 
DOUBLE  PRECISION  BLON( 1 ), BOUND, CANDA(COLMAX . PICK) . CANDCJC PICK) 
DOUBLE  PRECISION  CJ, COLA(COLMAX) , LJ, MEMR( LENR1 ) , UJ . UZBRO(MAIM) 
DOUBLE  PRECISION  ZBZERO(MAKM) , Y^(MAXM) , Z 
o 

o  Integer  arrays  and  variables 

o 

INTEGER  BNDTYP . COLLEN , DFEASQ , DTBRM , DUKBR , ERROR , FACTOR 
INTBGER  IOBRR , IOLOG . ITER , ITER 1 , ITER2 . LOOK . M , MAPI ( LENMI 1 ) 
INTEGER  MAPR( LENMR 1 ) , N , NTYPB2 , PRINT , TERMIN 
INTEGER  UNBDDQ , MEMI (LENIl).BASISC MAXM ) , CAND( PICK ) 

INTEGER  CANDlC COLMAX .PICK) , CANDLC PICK) , COLI ( COLMAX ) 

INTEGER  ROWTYP(MAXM) ,  STATUS (MAXN) 

O 

o  I/O  numbers  for  in  put,  error  output,  and  results  output 

o 

COMMON  /IQ_UNIT/IOIN ,  IOERR ,  IOLOG 
o 

o  Other  common  variables 

o 


o 


COMMON  /XMPLEN/  LENI,  LENMI.  LENMR.  LENR 
INTBGER  LENI,  LENMI,  LENMR,  LENR 


WA 


ss 


<nn 


Body  of  Program  (SUBROUTINE  YROBUST) 


LENI-LENI1 
LENMI-L2NMI1 
LENR-LENR1 
LBNMR-LENMR 1 
LOOK- 200 
FACTOR-50 
PRINT-0 
BNDTYP-1 

CALL  XMAPS ( BNDTYP , IOERR , MAPI , MAPR , MAXA , MAXM . MAXN . MEMI , MEMR ) 

Set  row  type  and  right  hand  sides  for  the  Z  equality 
constraints 

Nl-CR+l ) *NTYPE 
DO  100  1-1,  N1 
B(I)-0 
ROWTYP(I)-0 
CONTINUE 

Set  row  type  and  right  hand  sides  for  the  RO  and  ZST 
inequality  constraints 

N2-N1+A+A*NTYPB* (R+l ) 

DO  110  I-Nl+1 ,  N2 
B(I)-0 

RONTYP(I)-+l 

CONTINUE 

Set  row  type  and  right  hand  side  for  the  Y  equality  constraints 

N3-N2+NTYPB 
DO  120  I-N2+1,  N3 
B(I)-1 
ROWTYPC I )-0 
CONTINUE 

B ( N3+ 1 ) -REAL ( INT ) / TARGET  S 
R0WTYP(N3+ 1 )-0 

Let  N  he  the  current  number  of  variables;  incremented  by  XADDAJ 


Construct  the  first  NTYPE*S+1  structural  variables  associated 
with  the  defender's  strategy 


C-10 


DO  130  NT-1,  NTYPE 
DO  131  J-0,  S 
COLLBN-O 
DO  132  1-0,  R 

IP  (P(I,J)  .NB.  0)  THEN 
COLLEN-COLLEN+l 
COLA ( COLLEN ) - P ( I . J ) • V7RAC ( NT ) 

COLI (COLLEN ) -1+1+ (R+l ) * (NT-1 ) 

ENDIF 

132  CONTINUE 

COLLEN -COLLBN + 1 
COLA (COLLEN) -1.0 
COLI ( COLLEN ) -N2+NT 
IP  (J  .NB.  0)  THEN 
COLLEN-COLLEN+1 
COLA (COLLEN ) -  J • NFRAC ( NT ) 

COLI ( COLLEN )-N3+l 
ENDIF 
C 

c  Call  XADDAJ  to  enter  the  column  for  YJ  (or  variable  J+l) 

c 

OBJ-O 

CALL  XADDAJ ( OBJ . COLA , COLI , COLLBN . COLMAX . IOERR . 

1  X , MAPI , MAPR , MEMI , MEMR . N ) 

131  CONTINUE 

130  CONTINUE 

O 

c  Construct  the  columns  for  the  z  variables 

c 

DO  140  NT- 1 .NTYPE 
DO  141  I-l.R+1 
COLA(l)  — 1.0 
COLI( 1 )-I+(NT-l ) * (R+l ) 

DO  142  Il-l.A 

C0LA(I1+1)— 1.0 

COLKll  +  l  )-Nl+A+A*  (R+  1)*(NT-1)  +  (I1-1)*(R+1)  +  I 
142  CONTINUE 

c 

c  Call  XADDAJ  to  enter  the  column  for  z's 

c 

COLLEN- A+l 
OBJ-O 

CALL  XADDAJ ( OBJ , COLA , COLI . COLLEN , COLMAX , IOERR , 

1  X. MAPI, MAPR, MEMI. MEMR, N) 

141  CONTINUE 

140  CONTINUE 

C 

o  Construct  the  column  for  RO 

c 

DO  150  I-l.A 

COLA ( I ) - VBG ( I ) 

C0LI(I)-N1+I 
130  CONTINUE 


COLLEN-A 

OBJ-1 

CAliL  IADDAJCOBJ , COLA . COLI , COLLEN , COLKAX , IOERR , 

1  X , MAPI , HA PR , MEMI , MBMR , N ) 

C 

c  Construct  the  column  for  the  ‘s'  variables 

c 

DO  160  NT-1.  NTYPB 
DO  161  I-l.A 
COLA(l)  —  1 
COLIC 1)-N1+I 
DO  162  11-1,  R+l 
COLACIl+l )-l 

COLICII  +  D-NI+A+A'CR+D'CNT-D  +  CI-D'CR+D  +  II 
162  CONTINUE 

COLLEN -R 1-2 
OBJ-O 

CALL  IADDAJCOBJ , COLA . COLI , COLLEN , COLMAI . IOERR , 

1  X . MAPI . KAFR , MEN I , ME MR , N ) 

161  CONTINUE 

160  CONTINUE 

O 

c  Construct  the  column  for  the  ' t '  variables 

o 

DO  170  I-l.A 

COLAC 1 ) -REALCRVC I ) ) /REALC  TARGETS ) 

COLIC 1)-N1+I 
COLLEN- 1 

DO  171  NT-1 , NTT PE 
DO  172  11-1,  R 

COLLEN -COLLBN+ 1 


172 

171 

o 

o 

o 


170 

o 

o 

o 

ISO 

c 

o 

o 

o 

o 


COLA  C  COLLEN  )— 1 1  *  NFRAC  C  NT  ) 

COLICCOLXjBN)-N1+A+CNT-1) *A*  CR+1)  +  Cl-1) *  CR+1 )  +  Il*l 
CONTINUE 
CONTINUE 
OBJ-O 

Call  XADDAJ  to  enter  the  column  for  t  Cor  varl^le  S+3) 

CALL  XADDAJ  C  OBJ , COLA , COLI . COLLEN . COLKAX , IOERR . 

X . MAPI . KAPR , MEMI . ME MR . N) 

CONTINUE 

Start  all  the  structural  variables  at  their  lover  bounds 

DO  180  J-l,  N 
STATUSCJ)-0 

Set  the  number  of  constraints 
M-N3+1 

Construct  the  slack,  surplus  and  artificial  variables  with  XSLACX 


CALL  XSLACK(B , BASCB , BASIS . BASLB . BASUB , BLOW , 
BNDTYP , BOUND , 

COLA , COLI , COLMAX , IOERR , 

M ,  MAPI , MAPR , MAXM . HAXN , MEMI , MEMR . 

N ,  ROWTYP , STATUS , UZERO , XBZERO , Z ) 


Solve  the  LP  with  the  primal  simplex  method 


CALL  XPRIML(B , BASCB , BASIS . BASLB , BASUB , BNDTYP , BOUND , 
CAND . CANDA , CANDCJ . CANDI . CANDL , 

COLA , COLI . COLMAX . 

FACTOR , IOERR , IOLOG , ITER1 . I TER 2 . 

LOOK, 

M ,  MAPI , MAPR . MAXM , MAXN . MEMI . MEMR , 

N .  NTYPE2 , PICK , PRINT , STATUS . 

TERMIN , UNBDDQ , 

UZERO . XBZERO , YQ , Z ) 


Check  if  solution  has  been  found 


IF  (TERMIN  .NE.  1)  THEN 

WRITE  (IOERR, ' (A) ' )  ' 1 THERE  IS  AN  ERROR  IN  THE  FORMULATION 


WRITE  (IOERR, '(A, 14)') 
STOP 
ENDIF 


TERMINATION  CODE 


TERMIN 


Store  in  YII  the  solution  to  the  LP 


DO  200  I-l.NTYPE 
DO  201  J-0,  S 
YII(I,J)-0 
CONTINUE 
DO  210  I-  1,  M 

IF  (BASIS(I)  .LE.  (S+l ) *NTYPE)  THEN 

YII( (BASIS( I)-l ) / (S+l )+l ,MOD( BASIS ( I )-l,(S+l))) 

l-XBZERO(I) 

ENDIF 

CONTINUE 


RETURN 


Input  variable  type  declaration 

INTBGER  R,  S,  A,  1ST,  TARGETS 

INTEGER  MAT,  MAZE,  MAZS,  KAXNTYPE,  NTYPB 

INTEGER  RV(MAT) 

REAL  P(0 : MAZR , 0 : MAZS ) ,  YII(MAZNTYPE ,  0:MAZS),  NFRAC ( MAZNTYPE ) 
REAL  VFRAC( MAZNTYPE) 


ooo  o  o  o  oooooo  oo 


'"■in 


H 


« 


i 


S 


l 


c 
c 
G 
C 

C  VII(A) 

c 
c 
c 

c. •••••••• 

c 

c  Output  variable  type  declaration 

c 

REAL  IIKMAT.MAZNTYPE.O'.MAIR)  ,  VII(MAT) 

C 

C*  *  * . ••*••»*•••••• . . . ******** 

C 

c  Local  variables  needed  for  IMP.  See  IMP  Dictionary  for 

c  definitions, 

c 

c. . 

c 

c  Integer  parameters 

c 

INTEGER  MAUI. MAIN. MAIA.COLMAI. PICK 

PARAMETER  (MAIM-13.  MAIN-100.  MAXA-350.  COLMAI-3 .  PICK-6) 

INTEGER  LENI1 ,  LBNMI1 .  LENMR1 .  LENR1 
PARAMETER  (LBNI1-1000. LENMI 1 -9 , LENMR 1 -8 , LEHR 1-1000) 

Double  Precision  arrays  and  variables 

DOUBLE  PRECISION  B(MAXM) , BASCB(MAXM) . BASLB(MAIM) . BASUB(MAXM) 

DOUBLE  PRECISION  BLOW( 1 ). BOUND, CANDACCOLMAI, PICK) . CANDCJC PICK) 

DOUBLE  PRECISION  CJ, COLA(COLMAX) , LJ , MEMRCLENR1 ) , UJ . UZERO(MAIM) 

DOUBLE  PRECISION  XBZERO(MAXM) , YQ(MAZM) , Z 

Integer  arrays  and  variables 

INTEGER  BNDTYP , COLLBN . DFEASQ . DTERM . DUMBR , ERROR . PACTOR 
INTEGER  IOERR . IOLOG , ITER . ITER 1 , ITER2 , LOOK . M , MAPI ( LENMI 1 ) 

INTEGER  MAPR( LENMR 1 ) . N , NTYPE2 , PRINT . TERM IN 
INTEGER  UNBDDQ.MEMI(LENIl) .BASIS(MAXM) .CAND(PICK) 

INTBGBR  CANDKCOLMAX .PICK)  , CANDLC PICK)  , COLIC COLMAZ)  I 

INTEGER  ROWTYP(MAIM),  STATUS (MAIN) 

I/O  numbers  for  input,  error  output,  and  results  output 

COMMON/ I0_UNIT/ IOIN , IOERR , IOLOG 

* 

Other  common  variables  ^ 

COMMON /IMPLBN/LBNI,  LENMI,  LENMR,  LENR  ' 

INTBGBR  LENI,  LENMI.  LENMR,  LENR 


strategy  against  YII  for 
attaok  size  A:  fraction 
of  type  K  targets  assigned 
I  RV '  s 

REAL  The  expected  target  survival 

rate  associated  with  III 
and  YII 


< 


O  Q  t-  •-  t— *  ooooooo  *-»  o  o  o 


o 

c 

c 

c 

c 

c 


Body  of  Program  (SUBROUTINE  XROBUST) 


LENI-LBNI1 

LBNMI -LENMI 1 

LSNR-LENRl 

LENMR-LENMR1 

LOOK-BO 

PACTOR-BO 

PRINT-0 

BNDTYP-1 

CALL  XMAPSCBNDTYP . IOBRR . MAPI . MAPR . MAXA . MAIN . MAIN . HEM  I . MEMR ) 

Sat  rov  type  and  right  hand  side  for  the  equality  constraints 

DO  100  NT-1,  NTYPB 
B(NT)-1 
RO*TYP(NT)-0 
CONTINUE 

B ( NTYPE+ 1 ) -RBAL( RV(1 ) ) / REAL ( TARGETS ) 

Rorrrp(NTYPE+i)-o 

Let  N  be  the  ourrent  number  of  variables;  incremented  by  XADDAJ 
N-0 

Construot  the  NTTPB*(R+1)  structural  variables  associated  with 
the  attaoker's  strategy 

DO  110  NT-1,  NTYPB 
DO  111  1-0,  R 
COLLBN-1 
COLACD-l 
COLIC  1) -NT 
IP  (I  .HE.  0)  TB35N 
COLLEN-2 

COLA ( 2 ) - I • NPRAC ( NT ) 

COLI(2)-NTYPE+l 

KNDIP 

OBJ-O 

DO  112  J-0,  S 

OBJ--P( I , J ) * VPRAC(NT) * YII (NT , J)+OBJ 
CALL  XADDAJ ( OBJ , COLA , COLI , COLLSN , COLHAX , 

IOBRR . K . MAPI . MAPR . MEMI , ME MR , N ) 

CONTINUE 

CONTINUE 


Start  all  the  structural  variables  at  their  lover  bounds 


DO  12C  J-l  .  N 
STATUS ( J)-0 


Set  the  cumber  of  constraints 
M-NTYPE- 1 

Construct  the  slack  *>"rn’.us  and  artificial  variables  with  XSLACK 

CALL  IS LACK ( B . BASCB  BASIS  i \SLB . BASUB . BLOW , 

BNDTYP. BOUND 

COLA . COLI . COLMAI . IOBRR 

M .  MAPI . MAPR  MAXM  MAIN  MBMI . MBMR . 

N .  ROWTYP , STATUS  UZBRO . IBZERO . Z  ) 

Solve  the  LP  with  the  primal  simplex  method 

CALL  IPRIMLCB . BASCB, BASIS , 8ASLB . BASUB , BNDTYP . BOUND . 

CAND . CANDA . CAXDCJ , CANDI , CANDL . 

COLA . COLI . COLMAI , 

FACTOR . IOBRR , BO LOG . ITBR1 . ITBR2 . 

LOOK . M . MAPI . MAPR . MAXM . MAIN . MBMI . MBMR . 

N . NTYPB2 . PICK .PRINT , STATUS . TBRMIN , UNBDDQ , 

UZBRO, IBZERO. Y9.Z) 

Check  if  solution  has  been  found 

IF  (TBRMIN  NE.  1)  THEN 

WRITE  (IOBRR.  '  (A V)  ' 1 THERE  IS  AN  ERROR  IN  THE  FORMULATION 

WRITB  (IOBRR. '(A. 14)')  '  TERMINATION  CODE  -  TBRMIN 
STOP 
BNDIF 

Store  in  XII  and  VII  the  solutions  to  the  LP  for  RV(l) 

DO  200  I-  1,  NTYPE+I 

IF  (BASIS(I)  .LB.  NTYPB* (R+l))  THEN 
111(1.  ( BAS IS (I )-l)/(S+l)+l. 

MOD(BASIS(I ) -1 , ( S+ 1 ) ) ) -XBZERO( I ) 

ENDIP 

CONTINUE 

VII(l)—  Z 

Solve  for  the  remaining  RV's 

IF  (A  .Eg.  1)  THEN 
RETURN 
BNDIF 

DO  300  I  -  2.  A 

B  (  NTY PE -t  1 ) -RE AL (  RV  ( I )  )  /  RE AL (  T ARGET S  ) 


ooo 


1 

1 

1 


1 

1 

1 

1 

1 

1 


CALL  IBCOMP ( B . BASCB , BNDTYP , BOUND , 

COLA . COLI , COLMA1 , IOERR . 

M .  MAPI . MAPR , MAIM , MAIM , MEMI . MEMR . M , 

STATUS, IBZBRO.Z) 

CALL  IPRIMLCB , BASCB , BASIS , BASLB . BASUB . BNDTYP , BOUND , 
CARD , CANDA , CANDCJ , CANDI , CAMDL , 

COLA , COLI , COLMAI , 

FACTOR , IOERR , EOLOG . ITHR1 , ITER2 , 

LOOK , M . MAPI . MAPR , MAIM , MAIN . MEMI . MEMR . 

N .  NTYPE2 , PICI . PRINT . STATUS . TERMIN , UNBDDQ . 
UZERO , XBZERO , YQ , Z ) 


Store  in  III  and  VII  the  solutions  to  the  LP  for  RV(I) 


DO  301  II  -  1,  M 

IF  (BASIS(Il)  .LE.  NTYPE*(R+1))  THEN 
IIICI,  (BASIS(I1 )-l )/ (S+l )+l . 

1  MOD( BASISC II ) '1 . (S+ 1 ) ) )-IBZERO( II ) 

BNDIF 


301 

CONTINUE 

VII(I)—  Z 

300 

CONTINUB 

Q 

RETURN 


o 

END 


C-18 


.A-./  ■  r  . 


1 

1 


SUBROUTINE  SUMMARY ( NATTY PE , NDFTYPE , ANSWER , MINRV . MAXRV , 
INCRV , INT , TARGETS , R , S , PFA , PFD , PFD1 , 
PFD2 , OUT , MAXNTYPE , NTYPE , VTYPE , NTAR ) 


INTBGER  NATTYPE , NDFTYPE , MINRV , MAXRV , INCRV , INT , TARGETS , R , S 
INTEGER  OUT, MAXNTYPE, NTYPE 
INTEGER  NTARC MAXNTYPE) 

REAL  PFA , PFD , PFD1 . PFD2 , VTYPE (MAXNTYPE) 

CHARACTERS  ANSWBR . TITLE* 50 

WRITECOUT, ' (A) ' )  'OTHE  PARAMETERS  OF  THIS  PREALLOCATED  PREFERE 
1NTIAL  DEFENSE  GAME' 

WRITECOUT. 1000)  ‘O' 

TITLE- '0  THE  ATTACK  METHODOLOGY' 

IF  (NATTYPE  .BQ.  1)  THEN 

WRITECOUT, 1001)  TITLE,  SIMULTANEOUS' 

ELSE 

WRITECOUT, 1001)  TITLE,  ' SEQUENTIAL ' 

IF  (ANSWER  . BQ.  Y')  THEN 

WRITECOUT. 1001)  '  '  WITH  ATTACK  SIZE  AT' 

WRITBC OUT, 1001)  ’  ',  '  A  TARGBT  KNOWN  TO' 

WRITECOUT, 1001)  '  ',  '  THE  DEFENDER’ 

ELSE 

WRITECOUT, 1001)  '  '  WITH  ATTACK  SIZE  AT' 

WRITECOUT, 1001)  '  ' ,  '  A  TARGET  UNKNOWN  TO' 

WRITECOUT, 1001)  '  ',  '  THE  DEFENDBR' 

ENDIF 

ENDIF 

TITLE-'  THE  DEFENSE  METHODOLOGY' 

IF  (NDFTYPB  . EQ,  1)  THEN 

WRITECOUT, 1001)  TITLE.  'ONE  SHOT' 

ELSE 

WRITECOUT, 1001)  TITLE,  'SHOOT  LOOK  SHOOT' 

ENDIF 

TITLE- ' 0  THE  FAILURE  RATE  OF  RV  "  S ' 

WRITECOUT. 1002)  TITLE,  PFA 
IF  (NDPTYPE  .EQ.  1)  THEN 

TITLE-'  THE  FAILURE  RATE  OF  THE  INTERCEPTORS' 

WRITECOUT, 1002)  TITLE,  PFD 

ELSE 

TITLB- '  THE  FAILURE  RATE  OF  THE  FIRST  SALVO' 

WRITBCOUT , 1001 )  TITLE. 

TITLE-'  INTERCEPTORS' 

WRITBCOUT, 1002)  TITLE,  PFD1 

TITLE-'  THE  FAILURE  RATE  OF  THE  SECOND  SALVO' 

WRITBCOUT, 1001)  TITLE,  '  ' 

TITLE- '  INTERCEPTORS ' 

WRITECOUT, 1002)  TITLE,  PFD2 
BNDIF 

TITLB-' 0  MAXIMUM  NUMBER  OF  RV''S  ATTACKING  A  ' 

WRITECOUT. 1001)  TITLE,  '  ' 


TITLE- '  SINGLE  TARGET' 

WRITBC OUT, 1003)  TITLE.  R 

TITLE-'  MAXIMUM  NUMBER  OP  INTERCPTORS  DEFENDING  A 
WRITE (OUT, 1001)  TITLE,  '  ' 

TITLE-'  SINGLE  TARGET' 

WRITE(OUT , 1003)  TITLE,  S 

TITLE-’ 0  THE  MINIMUM  NUMBER  OF  RV''S' 

WRITECOUT, 1003)  TITLE.  MINRV 
TITLE-'  THE  MAXIMUM  NUMBER  OF  RV' 'S' 

WRITECOUT. 1003)  TITLE.  MAXRV 
TITLE- '  THE  ATTACK  SIZE  INCREMENT' 

WRITBCOUT, 1003)  TITLE.  INCRV 
TITLE-'  THE  NUMBER  OF  INTERCEPTORS' 

WRITECOUT. 1003)  TITLE.  INT 

TITLE-'  THE  TOTAL  NUMBBR  OF  TARGETS' 

WRITECOUT. 1003)  TITLE.  TARGETS 
TITLE- '0  THE  NUMBER  OF  TARGET  TYPES' 

WRITECOUT, 1003)  TITLE,  NTYPE 
DO  100  1-1,  NTYPE 

WRITECOUT. 1004)  '  TARGET  TYPE  '.  I.':' 

TITLE-'  NUMBBR  OF  TARGETS' 

WRITECOUT, 1003)  TITLE,  NTAR(I) 

TITLE-'  RELATIVE  VALUE' 

WRITECOUT, 1002)  TITLE.  VTYPB(I) 

CONTINUE 

WRITECOUT, 1000)  '0' 

FORMATC A , 72( ' • ' )) 

FORMAT (ABO,  A) 

FORMATC ABO,  FS.3) 

FORMATC ABO,  IB) 

FORMATCA, 12 . A) 


RETURN 

END 


SUBROUTINE  YRPRINTCYR,  MINRRV,  MAXRRV,  MAXS,  N,  OUT, 

1  MAINTYPE,  NTYPB) 

INTEGER  MINRRV.  MAXRRV.  MAXS.  MAXNTYPB ,  NTYPE 
REAL  YR( MAXNTYPB, 0: MAXS) 

INTEGER  OUT.  N1 ,  S 

WRITECOUT,  1000)  ' 1THE  ROBUST  DEFENSE  STRATEGY  FOR  RV  RANGE 
1.  MINRRV,  '  TO  ' .  MAXRRV, 

WRITECOUT, 1001)  '0' 

DO  100  NT-1, NTYPB 

WRITBCOUT, 1002)  ' 0 ’ . ' TARGBT  TYPE  ' , NT 
WRITBCOUT, 1008)  (I. 1-0,9) 

WRITECOUT, 1003)  ’  ' 


Nl-N/10 

IF  (HI  .  EQ.  0)  THEN 

WRITE(OUT, 1004)  (YR(NT, I) ,  1-0, N) 

ELSE 

WRITB(OUT, 1004)  (YR(NT , I) ,  1-0,9) 

DO  101  11-1,  N1 

IF  (II  ,EQ.  Nl)  THEN 

VRITE( OUT, 1004)  (YR(NT.I),  1-11*10,  N) 

ELSE 

WRITE(OUT, 1004)  (  YR  ( NT ,  I )  ,  1-11*10,11*10+9) 
ENDIF 
CONTINUE 
ENDIF 

WRITB( OUT, 1003)  '  ' 

CONTINUE 

WRITE (OUT. 1001)  '0' 

FORMAT( A,  15  ,  A ,  15  ,  A) 

FORMAT( A , 106( ' •  '  )) 

FORMAT(A, 471, A, 12) 

FORMAT(A, 71 , 93( ' - ' ) ) 

FORMAT ( 81. '  : ' , 10(F7 .4. '  :')) 

FORMAT(13K. 10(11,81)) 

RBTURN 

END 


SUBROUTINE  STPRINT ( STNAME , STBG , RV , A , N , OUT , MAZNAT . MAIN , 
1  MAZNTYPE , NTYPE , VFRAC , NFRAC ) 

INTEGER  MAZNAT. MAZN. MAZNTYPE, NTYPE 
CHARACTER *70  STNAME 

REAL  STBG (MAZNAT, MAZNTYPE. 0: MAZN),  NFRAC (MAZNTYPE) 

REAL  VFRAC (MAZNTYPE) 

INTEGER  RV(MAZNAT),  A.  N,  OUT 

WRITB(OUT. 1000)  STNAME 
WRITB(OUT, 1001)  '0' 

DO  100  NT-1,  NTYPE 

WRITB( OUT. 1000)  ’  ' 

WRITB(OUT, 1002)  'O',  ’TARGET  TYPE',  NT, 

1':  '.  100*NFRAC(NT) , '*  OF  TOTAL  TARGETS.  WITH  '.  VFRAC(NT) *  100 

l'%  OP  TOTAL  VALUB' 

IF  (NT  . BQ.  1)  THEN 

WRITB(OUT, 1007)  '  ATTACK  SIZE', (I,  1-0,9) 

BLSE 

WRITE(OUT, 1007)  '  ’,(1.  1-0,9) 

ENDIF 

WRITB( OUT. 1003)  '  ' 

DO  101  K  -  1,  A 
Nl-N/10 

IF  (Nl  . BQ.  0)  THEN 


ATTACK  SIZE', (I,  1-0,9) 
',(1.  1-0,9) 


WRITE(OUT, 1004)  RV(E),  (STBG(K,NT, I) .  1-0, N) 


102 


101 

100 


c 

1000 

1001 

1002 

1003 

1004 
1008 
1006 
1007 
C 


ELSE 

WRITE(OUT, 1004)  RV(K) .  (STBG(K, NT, I ) , 
DO  102  11-1.  N1 

IF  (II  .  EtJ .  Nl)  THEN 

WRITE(0UT, 1005)  (STBG(K,NT,  I) 

BLSB 

WRITE(OUT, 1005)  ( STBG(K , NT , I ) 

1  11*10+9) 

BNDIF 

CONTINUE 

ENDIF 

IF  (X  .NE.  A)  THEN 
WRITS(OUT, 1006) 


ELSE 

WRITE(OUT, 1003) 
ENDIF 
CONTINUE 
CONTINUE 

WRITB(OUT. 1000)  '  ’ 

WRITE (OUT, 1001)  ‘O' 


1-0,9) 


1-11*10,  N) 

1-11*10, 


FORMAT(A) 

FORMAT(A, 106( '• ' )) 

FORMAT ( A ,  17X ,  A12 , 2X ,  13  .  A .  F6 . 2 .  A ,  F6 . 2  .  A) 
FORMAT ( A. 1 06 ( '-')) 

FORMAT( 51 ,  15.  41,'  : ' . Z0(F7. 4 ,  ’  ‘  )) 

FORMAT ( 141, '  : ' , 10(F7 . 4 , '  :)) 

FORMAT ( 141 . 93 ( '-' )) 

FORMAT(A, 71, 10(11,61)) 


RETURN 

END 


SUBROUTINE  VPRINT( VNAME ,  VG,  RV.  A,  OUT.  MAXNAT) 
C 


o 

INTEGER  MAXNAT 
CHARACTER  *70  VNAME 
REAL  VG( MAXNAT) 

INTEGER  RV( MAXNAT),  A,  OUT 

Q 


WRITE (OUT, 1000)  VNAME 
WRITB( OUT, 1001)  'O' 

WRITE (OUT, 1000)  '0  ATTACH  SIZE 
WRITB(OUT, 1002)  '  ' 

DO  100  X-  1,  A 

WRITB(OUT, 1003)  RV(S),  VG(K) 
100  CONTINUE 

WRITE (OUT, 1002)  '  ' 


D-4 


r.n 


Utail* 


WRITE(OUT, 1001)  'O' 


o 

1000 

FORMAT(A) 

1001 

FORKAT(A,  29(  '  * 

')) 

1002 

FORMAT(A.29( 

')) 

1003 

FORMAT C  5X,  15, 

sx, 

c 

RETURN 

END 

:  ' .F11.4) 


SUBROUTINE  ALYRPRINT( YR ,  MINRRV ,  MAZRRV ,  MAZS,  N,  OUT, 
1  MAZNTYPE,  NTYPE ,  NFRAC . TARGETS ) 

c 

INTEGER  MINRRV.  MAZRRV,  MAZS,  MAZNTYPE ,  NTYPE 
REAL  YR ( MAZNTYPE, 0: MAZS) .NFRAC (MAZNTYPE) 

INTEGER  OUT.  Nl,  N,  TARGETS 

c 

WRITECOUT, 1000)  ' 1THE  ROBUST  DEFENSE  ALLOCATION  FOR  RV  RANGE 
1,  MINRRV,  ■  TO  ' .  MAZRRV,  ' 

WRITECOUT , 1001 )  '0' 

DO  100  NT-1, NTYPE 

WRITBCOUT, 1002)  ‘ 0 ', ' TARGET  TYPE  '.NT 
WRITBCOUT, 1005)  (I, 1-0,9) 

WRITECOUT, 1003)  '  ' 

Nl-N/10 

IP  (Nl  ,E Q.  0)  THEN 

WRITECOUT. 1004)  (YR (NT, I) 'NFRAC (NT) ‘TARGETS  ,  I-O.N) 

ELSE 

WRITECOUT. 1004)  (YR(NT, I) 'NFRAC(NT) 'TARGETS  ,  1-0,9) 
DO  101  11-1.  Nl 

IF  (II  .EQ.  Nl)  THEN 

WRITECOUT, 1004)  (YR(NT . I ) 'NFRAC(NT) 'TARGETS , 

1  I-I1'10,  N) 

ELSE 

WRITE  (  OUT ,  1 004  )  (  YR  (  NT .  I )  '  NFRAC  (  NT  )  *  TARGETS , 

1  I-I1'10, 11*10+9) 

ENDIF 

101  CONTINUE 

ENDIF 

WRITBCOUT, 1003)  '  ' 

100  CONTINUE 

WRITECOUT. 1001)  'O' 

o 

1000  FORMAT(A. IS , A , 13 , A) 

1001  FORMAT(A , 108(  "  ' ) ) 

1002  FORMATCA , 46Z, A , 12) 

1003  FORMATCA. 71, 93 ( )) 

1004  FORMAT ( 8Z, '  :',10(F7.1,'  :')) 

1003  FORMATC 13Z, 10(11 , 3Z) ) 

C 

RETURN 

END 


SUBROUTINE  ALPRINTC STNAME , STBG , RV , A , N , OUT , MAZNAT , MAIN . 
1  MAZNTYPE, NTYPE, VFRAC , NFRAC, TARGETS) 

c 

INTEGER  MAZNAT , MAIN , MAZNTYPE , NTYPE 
CHARACTER '70  STNAME 

REAL  STBG (MAIN AT, MAZNTY PE, 0: MAIN).  NFRAC (MAZNTYPE) 

REAL  VFRAC (MAZNTY PE) 


INTEGER  RV(MAXNAT),  A,  N,  OUT, TARGETS 

WRITECOUT, 1000)  STNAME 
WRITECOUT, 1001)  '0' 

DO  100  K-l.A 

WRITE (OUT, 1000)  '  ' 

WRITBC OUT, 1002)  'O'.  'ATTACK  SIZE  RV(K) 

IF  (K  .BQ.  1)  THEN 

WRITECOUT, 1007)  '  TARGET  TYPE',  (I,  1-0,9) 

ELSE 

WRITECOUT, 1007)  '  (I,  1-0  9) 

END  1 7 

WRITECOUT. 1003)  '  ' 

DO  101  NT-1, NTYPE 
Nl-N/10 

IF  (N1  .Eg.  0)  THEN 

WRITECOUT, 1004)  NT, 

1  ( STBGCK, NT, I )*NFRAC(NT) ‘TARGETS,  I-O.N) 

ELSE 

WRITBC OUT, 1004)  NT, 

1  (STBG(K,NT,I)*N7RAC(NT)*TARGETS,  1-0,9) 

DO  102  11-1,  N1 

IF  (II  . EQ.  Nl)  THEN 
WRITECOUT, 1003) 

1  ( STBG ( K . NT , I ) *  NFRAC ( NT ) *  TARGETS ,  1-11*10,  N) 

ELSE 

,  WRITECOUT, 1003) 

1  (STBG(K.NT.I)*NFRAC(NT)*TARGBTS,  1-11*10.11*10+9) 

BNDIF 

2  CONTINUE 

ENDIF 

IF  (NT  . NE .  NTYPB)  THEN 
WRITECOUT, 1006) 

ELSE 

WRITB (OUT, 1003)  '  ’ 

ENDIF 

1  CONTINUE 

0  CONTINUE 

WRITBCOUT, 1000)  ’  ' 

WRITB(OUT , 1001 )  '0' 

00  FORICAT(A) 

01  FORNATCA, 106( ' * ' )) 

02  FORNATCA. SOX. A13, 18) 

03  FORNATCA, 106C'-')) 

04  FORNATC8X,  13,  81,'  : ',10(P7.1.'  :')) 

08  FORNATC 141, *  10(77.1,'  :')) 

06  FORNATC 14X , 93( *  - ' ) ) 

FORNATC A. 7X, 10(11,81)) 


RETURN 

END 


SUBROUTINE  ALVPRINT( VNAME ,  VG,  RV,  A,  OUT,  MAXNAT, 
1  MAXNTYPE , NTYPE , TS ) 

c 
c 

INTEGER  MAXNAT , MAXNTYPE , NTYPE 
CHARACTER *70  VNAME,  CNAME,  BLANK 
REAL  VG( MAXNAT) ,  TS C MAXNAT , MAXNTYPE) 

INTEGER  RV( MAXNAT).  A,  OUT 
CHARACTER*10  DASH (20) ,  STAR(20) 
c 

BLANK-' 

WRITB(OUT, 1000)  VNAME 

c 

DO  200  NT-1,  NTYPE 

STAR(NT)- '•*•**••*«*' 

DASH(NT)-' - ' 

200  CONTINUE 

C 

CNAMB-BLANK( 1 : NTYPB*5 ) / / ' TARGET  TYPE' 

WRITE(OUT. 1001)  'O', (STAR(NT) ,  NT-1, NTYPE) 

WRITECOUT, 1002)  'O'.  CNAME 

WRITB (OUT. 1003)  ' ATTACK  SIZE ' ,  (I.  I-l.NTYPE) 

WRITE (OUT, 1004)  '  ' , (DASH(NT) ,  NT-1, NTYPB) 

DO  100  K-  1,  A 

WRITB(OUT, 1005)  RV(K) , (TS(K,NT) ,  NT-1, NTYPB) 

100  CONTINUE 

WRITE (OUT, 1004)  '  ' , (DASH(NT) ,  NT-1, NTYPE) 

WRITB(OUT, 1001)  'O', (STAR(NT) ,  NT-1, NTYPE) 

G 

1000  FORMAT(A) 

1001  FORMAT(A, 16( ' * ' ) , 10A) 

1002  FORMAT(A. 10X.A) 

1003  F0RMAT(2X.A11,2X, ' :  ',10(14,61)) 

1004  FORMAT ( A , 16 ('-'), 1 OA ) 

1005  FORMAT(5X,I5,5X, ' :  ' , 10(F7 . 2 , 3X) ) 

o 

RETURN 

END 
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SUBROUTINE  RVINTCOUNT( VNAME ,  TM,  RV.  A.  OUT, 

1  MAXNAT , MAXNTYPE , NTT PE , ROB ) 

o 
o 

INTEGER  MAXNAT  MAXNTYPE, NTYPE 
CHARACTER *70  VNAME ,  CNAME ,  BLANX 

INTBGBR  RV(MAXNAT),  A.  OUT,  TM ( MAXNAT , MAXNTYPE) .  N 
CHARACTER* 10  OASHC 20) ,  STAR(20) 

LOGICAL  ROB 

o 

BLANK-' 

WRITB(OUT, 1000)  VNAME 
DO  100  NT-1.  NTYPE 

STAR (NT)- '««•*•*••••' 

DASH(NT)- ' - ' 

100  CONTINUE 

CNAME-BLANXC 1 : NTYPE'3 ) / / 'TARGBT  TYPE' 

WRITB(0UT. 1001)  'O’, (STAR(NT) ,  NT-1 , NTYPB) 

WRITE (OUT, 1002)  'O',  CNAME 

WRITB( OUT, 1003)  ' ATTACK  SIZE ' ,  (I,  1-1, NTYPB) 

WRITE (OUT, 1004)  '  ' , (DASH(NT) ,  NT-1, NTYPB) 

IP  (ROB  ,BQ.  .FALSB.)  THEN 
DO  30  K-  1,  A 

WRITB(OUT, 1003)  RV(K) , (TM(K . NT) ,  NT-1, NTYPB) 


30 

CONTINUE 

BLSB 

WRITE(OUT, 1006)  'N/A',  (TM(1,NT) ,  NT-1, NTYPE) 

BNDIP 

WRITE (OUT, 1004)  '  ‘ 

1  , (DASH(NT) , 

NT-1, NTYPE) 

o 

1000 

1001 

1002 

WRITE (OUT. 1001)  '0 

PORMAT(A) 

FORMAT(A, 18( ' * ' ) , 10A) 
PORMAT(A.IOI.A) 

'  , (STAR(NT) , 

NT-1 .NTYPE) 

1003 

1004 

F0RMAT(2X,A11 ,2X, ' 
FORMAT(A, 18( ' - ' ) , 10A) 

10(14, ex)) 

1003 

FORMAT(3I,  15,  31.':  ' 

1 ,10(16,41)) 

1006 

O 

FORMAT (61,  A3,  6X.':  ' 

RETURN 

END 

1  , 10(16,41)) 
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PROGRAM  BATCH 


BATCH  creates  an  input  file  that  will  be  accepted  by  RPPDM 

PARAMETER  (MAXNTYPE-7 ,  MAXR-30,  MAES -30) 

INTEGER  R,  S , MAXRV , MINRV , INCRV , I NT , TARGETS 
INTEGER  NTAR ( MAXNTYPE ) 

REAL  VTYPE( MAXNTYPE) 


R 

THE 

S 

THB 

MAXRV 

THE 

MINRV 

THE 

INCRV 

THE 

INT 

THE 

TARGETS 

THE 

VTYPE 

THE 

NTAR 

THE 

MAXIMUM  NUMBER  OP  RV'S  AT  A  SINGLE  TARGET 
MAXIMUM  NUMBER  OF  INTERCEPTORS  AT  A  SINGLE 
TARGET 

MAXIMUM  ATTACH  SIZE 
MINIMUM  ATTACH  SIZE 
ATTACK  SIZE  INCREMENT 
NUMBER  OF  INTERCEPTORS 
NUMBER  OF  TARGETS 

RELATIVE  VALUES  OF  THE  TARGET  TYPES 
NUMBER  OF  TARGETS  FOR  EACH  TYPE 


INTEGER  NROBUST,  MAXRRV.  MINRRV 


NROBUST  THE  NUMBER  OF  ROBUST  DEFENSES  DESIRED 
MAXRRV  THE  MAXIMUM  ATTACK  SIZE  IN  THE  ROBUST  RANGE 

MINRRV  THE  MINIMUM  ATTACK  SIZE  IN  THE  ROBUST  RANGE 


INTEGER  TER . OUT 

CHARACTER* 12  FILE OUT,  NAME,  ANSWER* 1 


TER 

OUT 

NAME 

ANSWER 

FILBOUT 


THE  I/O  UNIT  NUMBER  FOR  THE  TERMINAL 
THB  I/O  UNIT  NUMBER  FOR  THE  OUTPUT  DEVICE 
A  CHARACTER  VARIABLE  SENT  BY  THE  USER 
A  CHARACTER  VARIABLE  SENT  BY  THE  USER 
NAMB  OF  THB  OUTPUT  DEVICE 


Set  TBR  to  default  I/O  unit  number  for  terminal 
TER-S 

Set  OUT  to  TER+1 
OUT-TER+ 1 

WRITECTER, ' (A) ' )  'OPlease  type  in  the  desired  file  name  (of  le 
lss  than  10  characters' 

WRITECTER, ' (A) ' )  '  Including  the  extension)  for  storage  of  t 
Ihe  parameters  ? ' 


ooo  ooo  ooo 


READ (TER, ' (A) ' )  FILEOUT 

OPEN(UNIT-OUT,  FILE-FILEOUT ,  STATUS- ' NEW ' ) 
c 
c 

c  Select  output  option 

c 

WRITE(TER , ' (A) ' )  ' lYou  have  three  options  for  the  output  of  th 
le  results  (with  RPPDM) :  ' 

WRXTECTER, '(A)')  '0  1)  TERMINAL  only' 

WRITE(TBR , ' (A) 1 )  '  2)  FILE  only' 

WRITB(TBR, ' (A) ' )  '  3)  TERMINAL  and  FILE' 

WRITE(TER , ' (A) ' )  ' OPlease  enter  the  number  for  the  desired  opt 
lion  ?' 

READ (TER,  FMT-* )  NOUT 
WRITE(OUT,FMT-* )  NOUT 

Ash  for  file  name,  if  necessary 

IF  (NOUT  . NE .  1)  THEN 

WRITE (TER, ' (A) ' )  'OPlease  type  in  the  desired  file  name  (o 
If  less  than  10  characters' 

WRITE(TER, ' (A) ' )  '  including  the  extension)  ?' 

READ (TER, ' (A) ' )  NAME 
WRITE(OUT.FMT-' (A12) ' )  NAME 
ENDIF 

Input  R  and  S 

WRITE (UNIT -TER,  FMT- ' (A . 12 , A) ' )  'OThe  MAXIMUM  number  of  RV '  '  s 
l(up  to  ' , MAXR, ' )  at  a  single  target  ?  ' 

READ(UNIT-TKR,  FMT-*)  R 
WRITE (UNIT-OUT,  FMT-*)  R 

WRITE (UNIT-TER,  FMT- ' ( A , 12 , A) ' )  ’OThe  MAXIMUM  number  of  INTERC 
1BPTORS  (up  to  '.MAXS,')  at  a  single  target  ?  ' 

READ (UNIT -TER,  FMT-*)  S 
WRITE (UNIT-OUT,  FMT-*)  S 

Select  attach  methodology 

WRITE(UNIT-TBR,  FMT-' (A)')  'OSeleot  one  of  the  following  attac 
lh  methodologies:  ' 

WRITE (UNIT-TER,  FMT- '(A)')  '0  1)  SIMULTANEOUS  ATTACK ' 

WRITE (UNIT-TER,  FMT-'(A)')  '  2)  SEQUENTIAL  ATTACK ' 

writb( UNIT-TER,  FMT- '(A)')  'OPlease  input  the  number  of  the  de 
1 sired  attaoh  ?  ' 

READ (UNIT-TER,  FMT-*)  NATTYPE 
WRITE (UN IT -OUT,  FMT-*)  NATTYPE 

c  Input  the  failure  rate  for  the  RV''s 

c 

WRITE (UNIT-TER,  PMT-'(A)')  'OThe  FAILURE  rate  of  the  RV"s  ?  ' 
RBAD( UNIT-TER,  FMT-*)  PFA 
WRITB( UNIT-OUT,  FMT-*)  PFA 
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Select  defense  methodology 

WRITE (UNIT-TBR,  FMT- '(A)')  'OSelect  one  of  the  following  defen 
lse  methodologies:  ' 

WRITE (TOIT-TBR,  FMT- '(A)')  '0  1)  ONE  SHOT' 

WRITE (UNIT-TBR ,  FMT- '(A)')  '  2)  SHOOT  LOOK  SHOOT' 

WRITB(0NIT-TBR,  FMT- '(A)')  'OPlease  input  the  number  for  the  d 
leslred  option  ?  ' 

RBAD( UNIT-TBR,  FMT-* )  NDFTYPE 
WRITB( UNIT-OUT,  FMT-* )  NDFTYPE 

If  the  the  attack  is  sequential,  find  out  whether  the  defender  knows, 
after  the  attaok  begins,  the  number  of  RV's  slated  for  each  target. 

IF  (NATTYPB  . SQ.  2)  THEN 

WRITE (UNIT-TER.  FMT- '(A)')  'OIs  the  defender  aware,  after 
lthe  attack  begins,  of  the  number ' 

WRITB(UNIT-TBR.  FMT-'(A)')  '  of  RV's  slated  for  each  targ 
let  (Y  or  N)  ?' 

RBAD( UNIT-TBR .  FMT- '(A)')  ANSWER 
WRITE (tJNIT -OUT ,  FMT-'(Al)')  ANSWER 
IF  (ANSWER  . BQ.  'Y')  THEN 

IF  (NDFTYPE  . EQ.  1)  THEN 

WRITB( UNIT-TER,  FKT- ' (A) ' ) ' ONOTE :  This  scenario  Is 
1  equivalent  to  one  vlth  a  simultaneous  attaok.’ 

BNDIF 

BNDIF 

BNDIF 

Find  the  failure  rates  for  the  interceptors 
IF  (NDFTYPE  . BQ.  1)  THBN 

WRITE (UNIT-TER.  FMT-'(A)')  'OThe  FAILURE  rate  of  the  Inter 
loeptore  ?  ' 

READ (UNIT -TER,  FMT-* )  PFD 
WRITE (UNIT-OUT.  FMT - * )  PFD 
BLSB 

WRITE(UNIT-TSR,  FMT-'(A)')  'OThe  FAILURE  rate  for  the  firs 
It  salvo  Interoeptors  ?  ' 

READ (UNIT-TER,  FKT-*)  PFD1 
WRITE (UNIT-OUT,  FMT-*)  PFD1 

WRITE ( UNIT-TBR ,  FKT- '(A)')  'OThe  FAILURE  rate  for  the  seco 
Ind  salvo  Interoeptors  ?  ' 

READ (UNIT-TER,  FMT-* )  PFD2 
WRITB (UNIT -OUT,  FMT-* )  PFD2 
BNDIF 

Specific  parameters  of  the  game 

WRITB(TBR, ' (A) ' )  'OThe  MINIMUM  attaok  size  ?' 

RE AD (TER,  FMT-*)  MINRV 

WRITE (UNIT-OUT.  FMT-*)  MINRV 

WRITB(TSR, ' (A) ' )  'OThe  MAXIMUM  attaok  size  ?' 

READ (TER,  FMT-*)  MAXRV 
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WRITBCUNIT-OUT,  FMT-* )  MAXRV 

WRITBCTBR, ' CA) ' )  'OThe  attack  size  INCREMENT  ?' 

READ(TBR,  PMT-* )  INCRV 
WRITBC UNIT-OUT,  FMT-* )  INCRV 

WRITBCTBR. ' (A) ' )  'OThe  NUMBER  of  Interceptors  ?' 

RBADCTBR.  fat-*)  INT 
WRITBC UNIT-OUT,  FMT-* )  INT 

WRITBCTBR, ' (A) ' )  ‘ OThe  TOTAL  NUMBER  of  targets  ?' 

RBADCTBR,  FMT-* )  TARGETS 
WRITBC UNIT-OUT.  PUT-*)  TARGETS 

WRITBCTBR, ' (A) ' )  'OThe  number  of  TYPES  of  targets  ?' 

RBADCTBR.  FMT-* )  NTYPE 
WRITBC UNIT -OUT,  FMT-* )  NTYPE 

C 

IF  CNTYPB  .NB.  1)  THBN 

WRITBCTBR,  'CA)')  ' OEnter  first  the  RELATIVE  VALUE  and  the 
In  the  NUMBER  of  targets ' 

WRITBCTBR.  'CA)')  ’  for  each  type.  Separate  the  two  entrie 
Is  for  each  target  type  with' 

WRITBCTBR,  ' CA) ' )  *  a  comma  and  hit  <CR>  following  the  ent 
lrles  for  each  target  type:  ' 
o 

o  Loop  through  eaoh  target  type 

o 

DO  100  I  -  1,  NTYPB 

RBADCTBR,  fmt-«)  VTYPBCI).  NTARCI) 

WRITBC UNIT-OUT,  FMT-* )  VTYPBCI).  NTARCI) 

TNT-TNT+NTARC I ) 

100  CONTINUE 

IF  CTNT  .NB.  TARGBTS)' THBN 

WRITBCTBR,  CA)')  'OThe  sum  of  the  targets  In  the  lndl 
lvldual  target  types  does  not ' 

WRITBCTBR,  'CA)')  '  equal  the  total  number  of  targets' 
STOP 
BNDIF 
END  IF 
O 
O 

WRITBCTBR. ' CA) ' )  lPlease  enter  the  number  of  different  ranges 
1  of  RV'  s  for  whloh  robust  ' 

WRITBCTBR, ' CA) ' )  ’  solutions  are  to  be  found  ?' 

WRITBCTBR. ' CA) ' )  '  •••••  Enter  0  if  no  robust  solution 

1  Is  de8l^ed*•••• 

RBADCTBR.*)  NROBUST 
WRITBC UNIT-OUT.  FMT-* )  NROBUST 
IF  C NROBUST  ,B9.  0)  THBN 
STOP 
ENDIF 

WRITBCTBR, ’CA)’)  OThe  lower  and  upper  bounds  for  the  RV  range 
Is  must  be  between' 

WRITBCTBR,  ' CIS  ,  '  1  and  M5)')  MINRV,  MAXRV 
o 

DO  200  I 1-1.  NROBUST 

WRITB  CTBR.'(A)')  '0  The  lower  bound  :' 
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READ  (TER,*)  MINRRV 

WRITE(UNIT-O0T,  FMT-* )  MINRRV 

NRITE  (TER. '(A)')  '0  The  upper  bound  :' 

READ  (TER,*)  MAXRRV 

WRITB(UNIT-OUT,  FMT-  * )  MAXRRV 

WRITE (TBR ,  '(//)') 

IP  (H  .NE.  1)  THEN 

WRITECTBR, '(A)')  '  NEXT1 
ENDIF 

WRITECTER, '(//)') 

200  CONTINUE 

C 

STOP 

END 


£ 


if 

P 


s', 


V 


I 


Let  K  be  the  set  of  all  target  types, 

T  be  the  total  number  of  targets, 

RV  be  the  number  of  reentry  vehicles  attacking  the  targets,  and 
INT  be  the  defending  interceptors. 

Let  R  be  the  maximum  number  of  RVs  allowed  to  attack  a  single  target  and 
S  be  the  maximum  number  of  interceptors  allowed  to  defend  a  single  target. 

LetVFk  =  ’VNk/  I_(Mc*Nk) 

keK 

and  NR  =  Nk/  I_  Nk 
keK 

where  Vk=  the  relative  value  of  a  type  k  target, 

Nk=  the  number  of  targets  of  type  k, 
and  k  is  a  target  type  in  K. 

Define  a  strategy  X  (or  Y)  as  the  set  of  (or  Y^)for  all  keK  where 

XM4 . xr) 

>Myo . r|) 


=  fraction  of  type  k  targets  assigned  i  RVs 
=  fraction  of  type  k  targets  assigned  j  interceptors. 


and 

4 


Let  P  be  the  matrix  [Pjj] 
and  F^the  ith  row  of  P, 

where  an  element  Pjj  is  the  probability  that  a  target  under  attack  by  i  RVs  and  defended 
by  j  interceptors  will  survive. 


The  formation  of  the  minimax  problem  is  as  follows. 


Let  VBG  =  expected  fraction  of  the  total  value  X_(Nk*  vk)  that  will  survive. 

keK 

VBG  =  max  ,™n  L  (VFk*Xk"  •  P’Yk) 

Yk>0  Xk>0  keK 


S  , 

I  Yf  =  I,  ke  K 

j=0  J 

s  v 

I  I  (j  •  Yf  •  NFk)  = 
keK  j=0 


INT 

T 


£  Xi  =1,  ke  K 

i=0 

keK  i=0  1 


where  X^  is  X^  transposed. 


G-1 


Taking  the  dual  of  the  inside  problem  yields  the  following  linear  program. 


LP1  is: 

VBG  =  max 

Y^O.sk.t 

subject  to 

Sk 

Sk  -  NFk- 1 


<  VT^*  Pq*  Yk 

<  vryPi*^ 


for  all  ke  K 


Sk-R»NFk;*t  ^VFjj'Pr'Y^ 


S  k 

X  Xj*YJ  *NFk=INT/T> 

keK  j-0 

Yj  >  0,  and  sk  and  t  unrestricted. 

The  Y  which  yields  VBG  is  the  optimal  minimax  defender's  strategy,  Y*. 


X*,  the  attacker's  minimax  strategy,  is  equivalent  to  the  dual  variables  of  the  in¬ 
equality  constraints. 


In  addition  s^  and  t  may  be  assumed  to  be  positive.  (See  [1]  on  Page  R-l.) 


2.  Robust  Game 

Assume  that  the  attacker  can  always  discover  and  thus  optimize  against  Y.  The 
defender  wishes  to  find  the  robust  strategy  YU  which  will  solve  the  following 
problem: 

max  min  (R^(X,Y)} 

Y=Yp  AeA 

where: 

(1)  Yp  is  the  set  of  all  Y's  such  that 

S  k 

g  Yj  =1  for  all  ke  K 

I_  I  G-nf^- Yk)  =-£!!. 
keK  J=° 
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A  ■-  V  -  A  ~  ^  ^  -V-  A AM  An  ^ 


(2)  A  is  the  set  of  possible  attack  sizes,  and 

(3)  Ra(X,Y)  =  VbG(A)  min 

the  set  ofallXk=(xo,...,XR) 


(V^.*Xk"*P  •  Yk) 


where 


£Xj  =  l  for  all  k  eK  and 
i=0 


X_  X  (i  *  NFk  *  Xf)  = 

teK  i=0 


k.  RV(A) 

i  /  —  rp 


Taking  the  dual, 
Ra(X,Y^=  VBG(A) 


subject  to 
sk(A) 


Yjf  >  0,  sk  (A),  t(A) 


^I_Sk(A)_RVIAIt(A)J 


^V^P0*^c 


sk  (A)  -  R  •  NF^'t  (A)  <  VIjj-Pi*  Yk 


for  all  keK.  . 


sk(A)  -  R  •  NFk* 1  (A)  -  VFk,pR*  ^ 


Inserting  this  back  into  the  original  problem  we  have: 

max  min  V  1  max  X  sk(A)~^f^  t  (A) 

YeYn  A  I  VBG(A)  sk(A),t(A)  keK 


which  yields  the  following  linear  program. 
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subject  to 


Y^>  0,  %(A),  t  (A) 


VBG(A)  p  <  X  sk(A)  — £V (A),  t  (A) 

keic  T 

sk(A)  -  VF^*  P()*Yk 

sk  (A)  -  NF  •  t  (A)  <  VI**  ?i  *\k 

sk(A)  -  R  •  NFfc*  t  (A)  ^  VPjr  •P^*'Vk; 


Once  more  we  can  assume  sj^A)  and  t  (A)  to  be  nonnegative  (see  [1]  on  Page  R-l). 
(Substituting  Z^(i)  for  VF**F}  *Yk  yields  the  LP  used  in  YROBUST.) 

By  assumption,  the  attacker  "knows"  YD,  the  Y  which  solves  LP2.  Thus  we  can  find 
XII,  the  optimal  attack  against  YTI,  and  VII,  the  expected  survival  rate  for  any 
attack  size  A,  by  solving  the  following  linear  program: 

LP3  is: 

VU(A)=  min  Z  (VFv«  Xk  "»P*  YIIk) 

Xk>  0  keK 

subject  to  p 

I  Xk  =1  for  all  kg  R 

r  i=0  1 

I  Z  i«Xjj  *NF,  =  RV(A) 
keK  i=0  T 


APPENDIX  H 

EQUATIONS  USED  IN  SIMAT1,  SIMAT2,  SEQAT1,  AND  SEQAT2 


TABLE  OF  CONTENTS 


1.  Equation  for  SIMATl-Simultaneous  Attack  with  One  Opportunity 

to  Shoot . 

2.  Equation  for  SIMAT2— Simultaneous  Attack  with  Two  Opportunities 

to  Shoot . 

3.  Equation  for  SEQAT1— Sequential  Attack  of  Unknown  Size  with  One 

or  Two  Opportunities  to  Shoot . 

4.  Equation  for  SEQAT2 --Sequential  Attack  of  Known  Size . 


Shoot 


I 

l 


We  assume  that  a  single  target  is  under  a  simultaneous  attack  by  A  identical 

missiles,  and  is  being  defended  by  D  identical  interceptors.  Let 

d  =  probability  that  a  defending  interceptor  will  destroy  the  attacking 
at  which  it  is  directed, 

a  =  probability  that  an  attacking  missile  will  destroy  the  target,  given  that 
it  evades  all  defending  interceptors. 

We  assume  that  the  defense  can  see  the  entire  attack,  and  must  decide  on  the  number  of 
interceptors  that  it  assigns  to  each  of  the  attacking  missiles. 

Given  suitable  independence  assumptions,  the  probability  of  the  target  surviving  an 
attack  of  nj  attacking  missiles,  each  of  which  are  being  attacked  by  j  defending  interceptors, 

is 

(1  -  a  (1  -d)))nj .  (1) 

Thus,  if  the  target  is  attacked  by  A  missiles  and,  for  each  j  between  0  and  D,  the  defender 
assigns  a  total  of  jnj  interceptors  against  nj  missiles,  such  that  j  interceptors  are  assigned 
against  each  of  these  nj  missiles,  where 

D 

I  nj  =  A 

j=0  J 

and 

D 

I  jn:=  D, 
j=0  J 

then  the  probability  that  the  target  survives  is 
D 

n  (1  -a(l -d)J)nj .  (2) 

j=o 


'Source:  Appendix  A  ot  Reference  [  1 1  on  page  R-l. 
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The  defender  wishes  to  select  the  nj  in  order  to  maximize  this  probability. 


We  wish  to  show  that  the  "uniform  defense"  obtained  by  spreading  the  D 
interceptors  as  equally  as  possible  over  the  A  attackers  is  optimal. 

Consider  an  allocation  of  interceptors  to  attackers  which  is  not  uniform.  Then  there 
is  a  pair  i  <  j  with  nj,  nj  >  0  where  i  +  2  £  j.  Consider  a  new  (and  more  uniform)  allocation 
obtained  by  allocating  i  +  1  interceptors  to  one  of  the  nj  attackers,  and  j  -  1  interceptors  to 
one  of  the  m  attackers.  The  probability  that  the  target  now  survives  is 

(1  -  a  (1  -  d)i+1)  (1  -  a  (1  -  d)M) 

(1  —  a  (1  -d)1)  (1  —  a  (1  —  d)J ) 

times  the  old  probability,  and  this  is  easily  shown  to  be  greater  than  1. 

The  most  uniform  of  defenses  assigns 
i  =  [D/ A]  defenders  to  n,  =  A(1-<D/A>)  attackers 

and 

j  =  [D/A]  +  1  defenders  to  nj  =  A  <D/A>  attackers 

(where  [x]  and  <x>  denote  the  integer  and  fractional  parts  of  x).  Thus,  the  optimal  defense 

is: 

{A(1-<D/A>)  forj  =  [D/A] 

A  <D/A>  forj  =  [D/A]  +  1 

0  otherwise. 

Substituting  these  values  for  nj  into  (2)  yields 

P(AJD)  =  (l-a(l-d)  [D/A]+i)A<D/A>(i_a(i_d)[D/A])A(I-<D/A>)> 

where  P(AJD)  is  the  probability  the  target  survives  given  that  it  is  attacked  by  A  missiles 
and  defended  by  D  interceptors  that  are  allocated  against  these  attacking  missiles  according 
to  this  uniform  defense.  Table  1  gives  the  numerical  values  of  P(A  J))  for  four  examples. 


Table  1.  VALUES  OF  P(A,D)  (SIMULTANEOUS  ATTACK) 
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The  following  material  is  extracted  from  Richard  M.  Soland,  "Optimal  Terminal 
Defense  Tactics  Against  Simultaneously  Arriving  RVs  When  Several  Sequential 
Engagements  are  Possible,"  TR-85/12,  Institute  for  Reliability  and  Risk  Analysis,  School 
of  Engineering  and  Applied  Science,  The  George  Washington  University,  18  October 
1985.  References  for  this  section  are  given  on  page  H-18. 

a.  INTRODUCTION 

An  important  and  frequent  element  of  ballistic  missile  defense 
models  is  the  presence  and  role  of  a  terminal  defense  that  defends  a 
single  target,  either  an  area  target  such  as  a  city  or  a  point  target 
such  as  an  ICBM  silo.  Given  that  such  a  terminal  defense  consists  of 
a  given  number  of  interceptors,  the  question  immediately  arises  as  to 
the  tactics  to  be  employed  by  the  defense  in  using  those  interceptors 
and  how  the  expected  damage  to  the  target  varies  as  a  function  of  the 
number  of  attacking  reentry  vehicles  (RVs). 

The  RVs  may  be  assumed  to  arrive  sequentially,  so  that  the  defense 
never  knows  how  many  are  coming,  or  else  simultaneously,  so  that  full 
knowledge  of  the  attack  size  is  obtained.  Intermediate  cases  may  be 
more  realistic  of  course.  The  defense  may  have  more  than  one  opportunity 
to  engage  some  or  all  of  the  RVs,  and  it  may  receive  information  on  the 
results  of  some  engagements  before  undertaking  subsequent  engagements. 
Here  we  treat  the  case  in  which  RVs  are  assumed  to  arrive  simultaneously 
and  a  fixed  number  of  engagements  of  each  RV  are  possible,  with 
shoot-look-shoot  capability • between  them. 

The  scenario  we  consider  is  as  follows:  a  single  target  is 


attacked  simultaneously  by  A  reentry  vehicles;  it  is  defended  by 
a  terminal  defense  consisting  of  D  interceptors.  The  RVs  act  independently, 

and  the  expected  fraction  of  the  target  destroyed  by  each  unintercepted  RV 
is  p,  where  0  <  p  <  1.  The  defense  may  engage  each  RV  up  to  K  times, 
with  one  or  more  interceptors  used  at  each  engagement,  and  it  has  shoot- 
look-shoot  capability  between  successive  engagements.  That  is,  the 
defense  observes  whether  or  not  a  particular  RV  survives  a  particular 
engagement  before  deciding  how  many  interceptors,  if  any,  to  use  against 
it  at  the  next  engagement.  All  interceptors  are  assumed  to  operate 
independently,  even  when  used  simultaneously  against  the  same  RV.  The 
single-shot  kill  probability  of  one  interceptor  against  one  RV  may  vary 
with  the  number  of  the  engagement  in  the  sequence  of  K  engagements. 

The  objective  of  the  defense  is  to  minimize  the  expected  fraction 
of  the  target  destroyed.  A  policy  for  the  defense  is  a  rule  that  specifies 
for  each  engagement,  how  many  interceptors  to  use  against  each  of  the  RVs 
that  has  not  been  destroyed  at  previous  engagements.  We  desire  to  find 
an  optimal  policy  as  a  function  of  A,  D,  and  K. 

Chapter  3  of  Eckler  and  Burr  (1972)  also  deals  with  the  tactics 
to  be  employed  by  a  terminal  defense  and  presents  several  different  models, 
each  based  on  specific  assumptions.  One  model,  that  of  section  3.5.1,  is 
a  special  case  of  ours;  it  corresponds  to  K  =  2  and  P  =  1.  Some  typical 
results  are  presented,  but  details  of  the  methodology  are  omitted. 

Burr  et  al.  (1985)  and  Falk  (1985b)  treat  the  case  in  which  RVs 
arrive  sequentially  but  only  one  engagement  of  each  RV  is  possible, 

'in ;  extends  such  analyses  to  the  case  of  several  engagements 


and  Fa  Ik  < I  i 


and  shoot-look-shoot  capability  between  them.  In  these  three  papers  the 
optimality  criterion  is  minimization  of  the  number  of  interceptors  D 
required  to  guarantee  that  the  expected  fraction  of  the  target  destroyed, 
as  a  function  of  A,  lies  below  a  given  bounding  function.  Our  criterion 
is  minimization  of  the  expected  fraction  of  the  target  destroyed;  as  we 
show,  however,  with  a  little  additional  work  we  can  also  treat  this  other 
criterion. 

A  brief  outline  of  the  paper  is  as  follows.  In  the  next  section 
we  define  notation  and  point  out  the  intuitive  result  (which  is  proved 
in  the  Appendix)  that  at  each  engagement  the  interceptors  used  should 
be  spread  as  uniformly  as  possible  among  the  attacking  RVs.  Using  this 
result,  we  present  a  dynamic  programming  algorithm  that  determines 
optimal  policies.  The  following  section  presents  several  extensions 
of  the  basic  model;  two  of  them  deal,  respectively,  with  the  expected 
number  of  interceptors  remaining  after  the  attack  and  determination  of 
the  minimum  number  of  interceptors  needed  to  provide  a  desired  level  of 
protection  of  the  target.  The  final  section  contains  illustrative 
numerical  examples. 

b.  ANALYSIS 

It  is  convenient  to  number  the  engagements  in  reverse  chronological 
order,  so  we  shall  refer  generally  to  there  being  k  engagements  remaining 
before  the  end  of  the  attack.  Here  k  =  0,1,..., K,  where  k  *  0  indicates 
that  no  further  engagements  are  possible.  For  k=l,...,K,  let  p^  be  the 
single-shot  kill  probability  of  one  interceptor  against  one  RV  on  the  kth 


engagement  before  the  end,  and  let  qR  =  1  -pfc.  We  assume  that  each  qk 
satisfies  0  <  qk  <  1;  the  contrary  cases  are  not  of  interest. 

We  define 

S(a,d,k)  *  the  expected  fraction  of  the  target  destroyed 
if  k  engagements  remain,  the  defense  has  d 
interceptors  left,  there  are  a  RVs  left 
undestroyed,  and  the  defense  follows  an 
optimal  policy  for  the  remaining  k  engagements. 

S(a,d,k)  is  defined  for  a=0,l,...,A;  d=0,l,...,D;  k=0,l,...,K.  Now 
define,  for  j=0,l,...,a;  a=0,l,...,A;  i=0,l,...,d;  d=0,l,...,D;  k-l,...,K, 

P(j|a,i,d,k)  =  the  probability  that  j  RVs  survive  the 

kth  engagement  from  the  end  when  there  are 
a  RVs  left  before  that  engagement  and  i 
of  the  d  remaining  interceptors  are  used 
in  an  optimal  manner  at  that  engagement. 

The  principle  of  optimality  of  dynamic  programming  now  allows  us  to 
write  the  following  recursive  equation  which  is  valid  for  a=l,...,A; 
d=0,l,...,D;  k=l,...,K: 

a 

S(a,d,k)  =  Min  -  {  Z  P(j|a,i,d,k)  S (j ,d-i ,k-l) } .  (1) 

i=0, . . . , d  j  =0 

In  order  to  use  this  recursion  to  calculate  S(a,d,k),  we  need  appropriate 
boundary  conditions  in  the  form  of  S(0,d,k)  and  S(a,d,0).  These  are 
clearly  as  follows: 


H-7 


S(0,d,k)  =  0  for  d=0,l . D;  k=l,...,K, 

S(a,d,0)  =  1  -  (1-p)  foi  a=0,l,...,A;  d=0,l,. 


.  .  ,D. 

It  remains  to  provide  the  P(j  ja,i,d,k)  before  the  above  recursion 
can  be  used.  This  is  rendered  relatively  straightforward  by  the  following 
result,  whose  proof  is  provided  in  the  Appendix:  an  optimal  manner  in 
which  to  use  i  of  d  interceptors  against  a  RVs  is  to  spread  them  as 
uniformly  as  possible  among  the  a  RVs.  Thus  if  i/a  is  an  integer,  say 
1,  each  of  the  a  RVs  is  assigned  I  interceptors  and  has  survival 
probability  .  The  number  of  RVs  that  survive  the  engagement  thus 
has  a  binomial  distribution,  so 

P(j|m,i,d,k)  -(*)  qj1  U-qJ)*'3  • 

More  generally,  i/a  is  not  an  integer.  Then  it  is  easily  shown  that 
a  +  a  LjL/aJ-i  of  the  a  RVs  are  assigned  [_i/aj  interceptors  each  and  the 
remaining  i-aLi/aJ  RVs  are  assigned  [~i/al  interceptors  each  (here  LxJ 
is  the  greatest  integer  less  than  or  equal  to  x  and  ["xl  is  the  smallest 
integer  greater  than  or  equal  to  x).  From  this  it  follows  that  the 
number  of  RVs  that  survive  the  engagement  is  the  sum  of  two  independent 
binomial  random  variables  with  slightly  different  success  probabilities. 
The  needed  probabilities  P(j|a,i,d,k)  may  then  be  obtained  by  numerically 
convoluting  the  two  binomial  distributions  in  question. 

An  alternative  way  to  obtain  the  needed  probabilities,  the  one 
actually  used  for  our  computational  results,  is  as  follows.  Let  J 
be  the  random  variable  who  probability  distribution  is  P (j ] a, i , d,k) . 


4 

j 

Then  J  is  the  sum  of  a  independent  Bernoulli  random  variables,  each  , 

I 

being  equal  to  one  if  a  specific  one  of  the  a  RVs  survives  the  engage-  ; 

ment,  and  equal  to  zero  otherwise.  We  may  thus  write  J  =  I£=1  X£,  ] 

where  P(X£  =  1)  =  r£  and  P(X£  =  0)  =  1  -  r£.  The  r£  are  given  by  j 

r  =  qa  for  2,=  l,...,a  +  aU/aJ  -  i  and  r£  =  q^  for  l  =  a  +  al_i/a]-i  +  l,...,a, 

where  a  =  Li/aJ  and  6  =  fi/al .  By  defining  the  random  variables 

j  E  £S  x  for  s  =  1, . . . ,a,  so  that  J  =  J  ,  we  can  calculate  the 
s  £=i  a  a 

probability  distributions  of  the  Jg  successively  from  the  recursion 
(which  is  valid  for  j=l,...,s  and  s=2,...,a) 

P(Js  -  J)  -  U-rs)  PCJ^M)  + 

along  with  the  boundary  conditions 

s 

P(J  «  0)  -  n  (l-rc) ,  P(J  -1)  =  r  . 

S  A=1  *  1  i 


A  closed-form  expression  for  S(a,d,l)  is  sometimes  useful,  and 

it  is  easily  obtained  from  the  uniform-defense  property  and  fairly 

simple  probabilistic  analysis;  the  result  is 

c/  J  !  fi  „  a\ a+aa-d  „  SNd-aa 

S (a ,d , 1)  =  1  -  (1  -  pq^)  (1-pq  )  , 

where  a  =  Ld/aJ  and  8  =  fd/al. 


C-  EXTENSIONS  OF  THE  BASIC  MODEL 

In  anticipation  of  the  possibility  of  another  attack  after  the 
present  one  is  over,  the  terminal  defense  may  be  Interested  in,  as  a 
secondary  criterion,  the  expected  number  of  interceptors  remaining  to 
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it  at  the  end  of  the  current  attack.  We  thus  define 

T(a,d,k)  =  the  expected  number  of  interceptors  left  after 
the  attack  if  k  engagements  remain,  the  defense 
presently  has  d  interceptors  left,  there  are 
a  RVs  left  undestroyed,  and  the  defense  follows 
an  optimal  policy  for  the  remaining  k  engagements. 

Like  S (a ,d ,k) ,  T(a,d,k)  is  defined  for  a=0,l, — ,A;  d=0,l,...,D; 
k=0, 1 , . . . , K.  We  can  calculate  T(a,d,k)  along  with  S(a,d,k)  as  follows: 

Let  i*  be  a  minimizing  value  of  i  in  the  recursion  (1)  used  to  determine 
S(a,d,k) .  Then 

a 

T(a,d,k)  =  l  P(j  |a,i*,d,k)  T(j ,d-i*,k-l) ,  (2) 

j=0 

for  a=l,...,A;  d=0,l,...,D;  k=l,...,K.  The  necessary  boundary  conditions 
are 

T(0,d,k)  =  d  for  d=0,l,...,D;  k=0,l,...,K, 

T(a,d,0)  =  d  for  a=l,...,A;  d=0,l,...,D. 

The  recursion  (2)  does  not  serve  to  uniquely  determine  the  T(a,d,k) 
unless  the  minimizing  i*  in  (1)  are  unique.  But  the  index  i*  is  not 
necessarily  unique,  so  we  adopt  the  following  convention,  which  then 
serves  to  determine  a  unique  value  for  T(a,d,k).  If  i*  in  (1)  is  not 
unique,  use  the  smallest  value  of  i*  that  maximizes  T(a,d,k)  as  computed 
by  (2). 

The  boundary  condition  S(a,d,0)  =  1  -  (l-p)a  given  above  was  based 
on  the  assumption  of  RVs  that  act  independently  and,  if  unintercepted. 


each  destroy  an  expected  fraction  p  of  the  target.  The  dynamic 
programming  scheme  given  works  equally  well  with  the  more  general 
boundary  condition  S(a,d,0)  =  g(a)  for  a=Q,l,...,A;  d=0,l,...,D,  where 
g  is  nondecreasing,  g(0)  =  0,  and  g(A)  <  1.  This  allows  use  of  an 
arbitrary  relationship  between  the  number  of  unintercepted  RVs  and  the 
expected  fraction  of  the  target  destroyed. 

In  damage-limitation  studies  it  is  sometimes  of  interest  to 
determine  the  smallest  number  of  interceptors  needed  by  the  defense  to 
provide  a  desired  level  of  protection  of  the  target.  One  way  to 

interpret  the  phrase  "desired  level  of  protection"  is  by  specifying  a 
nondecreasing  maximum  damage  function  f  and  requiring  that  the  expected 
fraction  of  the  target  destroyed  by  a  RVs  not  exceed  f(a)  for  a=l,...,A. 

For  example,  see  Burr  et  al.  (1985)  and  Falk  (1985b)  for  analyses  of  this 
problem  when  the  RVs  are  assumed  to  arrive  sequentially  and  the  defense 
has  no  shoot-look-shoot  capability,  and  see  Falk  (1985a)  for  an  extension 
to  the  case  in  which  the  defense  does  have  shoot— look— shoot  capability. 

In  the  present  context,  and  for  fixed  A  and  K,  we  may  phrase  the  problem  as 

Minimize  D 

subject  to  S(a,D,K)  <  f(a),  a=l,...,A. 

We  can  solve  this  problem  easily  by  continuing  the  computations  of  the 
dynamic  programming  scheme  for  successively  larger  values  of  D  until 
one  is  found  that  satisfies  all  the  indicated  constraints. 


d.  EXAMPLES 


For  the  case  p  =  0.7,  =  0.8  and  =  0.9,  Table  1  gives 

S (A, D, 2) ,  i*,  and  T(A,D,2)  for  A,D  =  1(1)10;  the  three  quantities 
appear  in  respective  rows.  For  example,  with  A  =  7  RVs,  D  =  9  interceptors 
and  2  engagements  left,  the  defense  should  use  i*  ■  7  interceptors  and 
hold  2  in  reserve  for  possible  use  at  the  last  engagement.  The  expected 
fraction  of  the  target  destroyed  is  S(7,9,2)  =  0.0631  and  the  expected 
number  of  interceptors  left  after  the  attack  is  T(7,9,2)  =  0.957. 

Figure  1  shows  the  expected  fractional  damage  S(A,D,2)  as  a 
function  of  A  for  D  =  2(2)10. 

Figure  2  shows  solutions  of  the  damage  limitation  problem 
minimize  D  subject  to  S(a,D,K)  <  f (a) ,  a“l,...,10,  for  the  linear 

function  f(a)  =  sa.  Solutions  are  given  for  a  range  of  values  of  the 
slope  s  and  for  both  K  ■»  1  and  K  *■  2;  the  other  parameters  of  the 
problem  are  P  ■  0.9  and  p^  *  *  0.7.  For  example,  if  it  is  desired 

to  have  a  defense  whose  expected  fractional  damage  per  attacking  RV 
is  limited  to  s  3  0.125,  then  9  interceptors  are  needed  if  only  one 
engagement  is  possible  (K  **  1),  but  only  8  interceptors  are  required 
if  two  engagements  are  possible  (K  ■  2) . 
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TABLE  1 


Results  for  2  Engagements  when 
*  0.8,  p£  *  0.9  and  p  -  0.7 


TABLE  OF  EXPECTED  FRACTION  0T  TARGET  DESTROYED,  OPTIMAL  NUMBER  OF 
INTERCEPTORS  TO  USE.  AND  EXPECTED  NUMBER  OF  INTERCEPTORS  LEFT  OVER 


0  =  1  0.0700  0.7210  0.9163  0.9749  0.9925  0.9977  0.9993  0.9998  0.9999 
D=  1  1  1  1  1  1  1  1  1  1 
D=  1  0.000  0.000  0.000  0.000  0.000  0.000  0.000  0.000  0.000 


1  0  =  A 

1 . 0000 
1 

0.000 


0  =  2  0.0070  0.1351  0.7405  0.9222  0.9766  0.9930  0.9979  0.9994  0.9998  0.9999 
0=  2  2  22  2  2  2  22  2  2 
0  =  2  0.000  0.000  0.000  0.000  0.000  0.000  0.000  0.000  0.000  0.000 

D  =  3  0.0007  0.0326  0.1  956  0.7537  0.3276  0.9733  0.9935  0.99S0  0.9994  0.9996 
0=  3  3  2  3  3  3  3  3  3  3  3 
D=  3  0.000  0.810  0.000  0.000  0.000  0.000  0.000  0.000  0.000  0.000 

D=  4  0.0001  0.0076  0.0550  0.2519  0.7756  0.9327  0.9798  0.9939  0.9962  0.9995 
0=  4  4  2  3  4  4  4  4  4  4  4 
D=  4  0.000  1.620  0.729  0.000  0.000  0.000  0.000  0.000  0.000  0.000 

D=  5  0.0000  0.0026  0.0146  0.0603  0.3043  0.7913  0.9374  0.9812  0.994d  0.Q96? 
0=  55234555555 
D=  5  0.000  2.430  1.453  0.656  0.000  0.000  0.000  0.000  0.000  0.000 

D=  6  0.0000  0.0006  0.0062  0.0237  0.1079  0.3530  0.8059  0.9418  0.9825  0.9943 
0=  6  6  4  3  4  5  6  6  6  6  6 
0=  6  0.000  1.960  2.1S7  1.312  0.590  0.000  0.000  0.000  0.000  0.000 

D=  7  0.0000  0.0001  0.0020  0.0110  0.0349  0.1373  0.3983  0.6195  0.9458  0.9638 
D=  7  7  4  3  4  56  7  7  7  7 
D=  7  0.000  2.940  2.916  1.968  1.181  0.531  0.000  0.000  0.000  0.000 

D  =  8  0.0000  O.OOCO  0.0009  0.004  1  0.0171  0.043  1  0.1  679  0.4404  0.8321  0.9496 
D  =  6  8  4  4  4  5  6  7  8  8  6 
D=  8  0.000  3.920  3.208  2.624  1.771  1.063  0.478  0.000  0.000  O.OOC 

0=  9  0.0000  0.0000  0.0002  0.0024  0.0069  0.0245  0.0631  0.1994  0.4796  0.8439 
D  =  9  9  4  6  4  5  6  7  3  9  9 
D=  9  0.000  4.901  2.911  3.260  2.36?  1.594  0.957  0.430  0.000  0.000 

D  = 1 0  0.0000  0.0000  0.0000  0.0009  0.004?  0.0105  0.0330  0.0799  0.231  3  0.5160 
0=10  10  4  6  4  5  6  7  8  9  10 
0=10  0.000  5.331  3.881  3.937  2.95?  2.126  1.435  0.861  0.367  0.000 
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ANNEX 


Here  we  provide  a  verification  of  the  intuitive  result  that, 
regardless  of  the  number  of  engagements  remaining,  an  optimal  manner 
in  which  to  use  i  of  d  interceptors  against  a  RVs  is  to  spread  them  as 
uniformly  as  possible  among  the  a  RVs,  i.e.,  use  a  uniform  defense. 

Suppose  there  are  k  >  1  engagements  remaining  and  suppose  there 
is  an  optimal  way  to  use  i  of  d  interceptors  against  the  a  RVs  that  is 
not  uniform,  i.e.,  there  are  two  RVs  such  that  one  of  them  is  assigned 
at  least  2  more  interceptors  than  the  other.  We  will  show  that  the 
expected  damage  done  to  the  target  is  not  increased  if  one  interceptor 
is  switched  from  the  heavily  attacked  RV  to  the  lightly  attacked  RV. 

By  repetition  of  this  step  a  finite  number  of  times  it  then  follows 
that  a  uniform  defense  is  optimal. 

Let  u(£)  be  the  number  of  interceptors  assigned  to  RV  l  (£.  =  1 , . .  . , a)  , 
so  that  Z ^  u(£)  =  i.  Without  loss  of  generality,  we  may  assume 
u(l)  >  u(2)  +  2.  Define  random  variables  X^,...,Xa  as  =  1  if  RV  Z 
survives  its  attack  by  u(£-)  interceptors  and  *  0  otherwise. 

Let  X, _  =  X,  +  X„.  Then  J  ;  X.,  +  .  X  is  the  number  of  RVs 

12  1  2  12  £=3  l 

that  survive  the  current  engagement  and  E[S(J,  d-i,  k-1)]  is  the  expected 
fractional  damage  done  to  the  target. 


Now  consider  the  alternative  defense  defined  by  u  (1)  =  u(l)  -1, 

u'(2)  3  u(2)  +  1,  and  u’(£)  =  u(£)  for  l  >  2.  Define  the  random  variables 

X',...,X*  in  the  same  manner  as  before,  and  let  X'  =  X'  +  X'.  Then 
la  12  1  2 


J'  =  X*  +  Zn  .  X'  and  E(S(J',  d-i,  k-1)]  are  interpreted  as  above,  but 


for  the  alternative  defense  instead.  X'  ,  X',...,X'  are  mutually 

J  £1 

independent,  as  are  X,  »,  X-,...,X  . 

iz  J  a 

st 

We  now  show  that  X|2  <  X^2,  ^'e,»  xj_2  stochasti cally  less 

than  X^?  [see  Barlow  and  Proschan  (1975),  p.  110].  X^  and  X^2  have 
possible  values  0,  1  and  2,  and 


P(X12=0)  =  (l-q£(1))  (l-q^(2))  =  1  -  qJJ(1) 


u(2)  u(l)+u(2) 

qk  +  qk 


P(x12=i) 


,  u(l)+u(2) 

2  qk 


P(X 


u(l)+u(2) 

qk 


Corresponding  expressions  hold  for  the  probability  distribution  of  X^2> 
with  u’ (1)  =  u(l)-l  and  u’ (2)  =  u(2)  +  1  substituted.  Since 
u(l)  +  u ( 2)  =  u’(l)  +  u ' (2) ,  it  follows  that  P(X  >  1)  =  P(X^2  >  1). 
Simple  algebra  yields 


0)  -  PCX. 


u(2) 

qk 


1)  (1 


u ( 1 ) -u (2)-l. 
9l.  } 


and  this  suffices  to  show  that  X^2  S£  X^.  Since  x£  =  X£  f°r  ^  >  2, 

it  follows  from  a  simple,  extension  of  exercise  1  on  page  176  of  Barlow 

st 

and  Proschan  (1975)  that  J'  <  J.  Since  S(j ,  d-i,  k-1) ,  by  virtue  of 
its  definition,  is  a  nondecreasing  function  of  j  for  fixed  values  of 
the  other  arguments,  it  follows  from  exercise  15  on  page  151  of  Barlow 
and  Proschan  (1975)  that  E[S(J\  d-i,  k-1)]  <  E[S(J,  d-i,  k-1 ) j .  Thus 
the  alternative  defense  does  not  increase  the  expected  damage  done  to 
the  target,  and  our  verification  is  now  complete. 


REFERENCES 


[1]  BARLOW,  R.  E.  and  F.  PROSCHAN  (1975).  Statistical  Theory  of 

Reliability  and  Life  Testing.  Holt,  Reinhart  and 
Winston,  New  York. 

[2]  BURR,  S.  A.,  J.  E.  FALK  and  A.  F.  KARR  (1985).  Integer  Prim-Read 

Solutions  to  a  Class  of  Target  Defense  Problems.  Operations 
Research ,  33,  726-745. 

[3]  ECKLER,  A.  R.  and  S.  A.  BURR  (1972).  Mathematical  Models  of 

Target  Coverage  and  Missile  Allocation.  Military 
Operations  Research  Society,  Alexandria,  VA. 

[4]  FALK,  J.  E.  (1985a).  Prim-Read  Defense  with  Defender  Shoot-Look- 

Shoot.  Memorandum.  Institute  for  Defense  Analyses, 
Alexandria,  VA. 

[5]  FALK,  J.  E.  (1985b;.  Minimal  Defenses  Enforcing  Arbitrary 

Maximum  Damage  Bounds.  Proceedings  of  the  1985  IEEE 
International  Conference  on  Systems ,  Man  and  Cybernetics . 
Tucson,  Arizona,  November  12-15,  1985,  651-655. 


I 

r 


J 

4 


I 

» 


There  is  a  single  target  that  may  come  under  an  attack  by 
an  unknown  number  of  sequentially  arriving  RV's. 

We  are  protecting  this  target  with  D  interceptors,  and  we 
have  sufficient  time  to  observe  the  results  of  a  first  volley 
against  an  RV,  and  fire  a  second  volley,  if  necessary. 

Let 

q  =  probability  that  a  defender  targeted  against  an  RV 
will  miss  in  the  first  volley, 

r  =  probability  that  a  defender  targeted  against  a  surviving 
RV  will  miss  in  the  second  volley, 

s  =  probability  that  an  RV  surviving  both  volleys  will  fail 
to  destroy  the  target,  and 

p(A,D)  *  probability  that  the  target  protected  by  D  defenders 

Is  destroyed  by  A  weapons,  given  some  firing  doctrine. 

Figure  1  illustrates  a  pair  of  possibilities  when  D=l. 

Here  (q,r,s)  =  (0.1,  0.2,  0.3),  and  the  results  of  3  firing; 
doctrines  are  displayed.  The  doctrines  are 

a)  Don't  fire  at  any  RV  (unprotected  case  =  D=G), 

b)  Fire  one  defender  in  the  first  volley  against  the 
first  arriving  RV. 

c)  Fire  one  defender  in  the  second  volley  (none  in  the 
first)  at  the  first  arriving  RV. 


“‘‘Source:  James  -  .  Falk,  "Prim-Read  Solutions  with  Shoot-Look- 

Shoot",  uncublished  memorandum,  Institute  for  Defense  Analyses, 
31  July  1535. 


H-19 


S'l 


Figure  1.  Expected  Damage  Curves  for 
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With  only  one  defender,  and  with  its  chances  of  success  being 
higher  if  used  in  the  first  volley,  it  is  clear  that  it  should 
be  so  used.  With  several  defenders,  and  different  q's  and  r’s, 
the  situation  becomes  less  clear. 

In  keeping  with  the  Prim-Read  philosophy,  we  will  use  the 
maximum  of  the  ratios  p(A,D)/A  ( A=1 , 2 , 3 , • . • )  as  the  measure  of 
effectiveness  of  a  firing  doctrine.  In  Figure  1,  we  have 

maximum  ratio  with  doctrine  (a)  =  0.7  (p  (1,1)/1) 

maximum  ratio  with  doctrine  (b)  =  0.361  (p^(2,l)/2) 

maximum  ratio  with  doctrine  (c)  =  0.371  (p^(2,l)/2) 

so  that  firing  doctrine  (b)  is  the  most  effective  in  that  its 
maximum  ratio  is  the  smallest. 

In  this  note,  we  will  assume  that  the  defender  will  always 
decide  on  his  first/second  volley  allocations  in  such  a  way  as 
to  minimize  the  maximum  of  those  ratios,  where  the  minimization 
takes  place  over  all  possible  such  allocations.  Under  this 
"behavioral  assumption,”  we  make  the  following  definitions. 

For  each  D=0,l,2,...,  let 

Q, (A,D)  =  probability  that  the  target  is  destroyed  by  A 

attackers  given  that  there  are  D  defenders  and  k 
of  them  are  sent  against  the  first  arriving  RV 
in  a  first  volley,  with  a  second  possible,  if 
needed , 

R^(A,D)  =  probability  that  the  target  is  destroyed  by  A 

attackers  given  that  there  are  D  defenders,  the 
defense  has  but  a  single  (second)  volley  at  the 
first  arriving  RV  and  sends  k  defenders  at  it, 


when,  in  each  case,  k  ranges  from  0  to  D. 


For  any  given  D,  set 


m(D)  =  min  max  Qk(A,D)/A  (I) 

0<k<D  Ae I+ 

and  let  k*(D)  denote  the  smallest  Integer  such  that  the  above 
minimum  is  obtained.  Given  the  functions  Q^(*,D),  Q^(*,D),..., 
3n(*,D),  both  m(D)  and  k*(D)  are  well-defined,  and  k*(D) 
represents  the  number  of  interceptors  that  the  defender  fires 
in  his  first  volley  against  the  first  RV  that  he  sees  when  he 
has  D  interceptors  left. 

In  the  event  that  a  defender's  first  volley  fails,  he  is 
left  with  a  "reserve  volley".  Assuming  that  he  now  has  D 
defenders  left,  if  we  have  the  functions  RQ(*,D)  R^(*,D),..., 

Rp( • , D) ,  we  define 

n(D)  =  min  max  R^(A,D)/A  (2) 

0<k<D  Ae I+ 

and  let  k**(D)  denote  the  smallest  integer  which  minimizes  the 
above.  Then  k**(D)  represents  the  numbers  of  interceptors  that 
the  defender  would  fire  if  he  had  but  one  (second)  volley  to  shoot 
at  the  first  RV  that  he  sees  with  D  interceptors  left  (he  has 
two  volleys  at  any  subsequent  RV),  and  he  wishes  to  choose  a 
damage  curve  with  the  smallest  slope. 

Let 


(3) 

( a ) 

For  each  D,  these  represent  the  actual  expected  damage  curves 
that  the  defender  has  selected. 


Q(A,D)  =  Qk*(D)(A,D) 

and 

R(A,D)  -  R^#  #  ( r; )  (  A ,  D  ) 
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Note  that 


Q  (  A ,  0  )  =  R  ( A  ,  0  )  =  l-sH  A=0 , 1 , .  .  .  . 


(if  the  attackers  are  perfect,  this  holds  if  we  define  0J=1). 

The  following  recursions  hold 

R,  ( A ,  D )  =  rk ( 1-s )  +  ( l-rk ( 1-s ) )  Q(A-l,D-k)  (5) 

Qk(A,D)  =  qkR ( A , D-k )  +  (l-qk)  Q(A-l,D-k)  (6) 

for  k=Q,l,...,D  and  D^_l .  They  may  be  solved  in  the  order: 

Q ( A , 0 )  =  R( A, 0 )  -  R1(A,1)  =  R(A,1) 

R(A,1)  -  Qq(A,1),  Q1(A,1)  -  Q(A, 1) 

Q( A , 1 )  -  R1(A,2),  R2(A,2)  -  R( A, 2 ) 

R(A,2)  -  Q3(A,2),  Q1(A,2),  Q2(A ,2)  -  Q(A,2)  etc. 

Note  that  R0(A,D)  is  never  computed.  The  above  recursion  for 
Rq(A,D)  is 

"Rq ( A , D )  =  (1-s)  +  s  Q(A-1,D) 

and  the  function  Q(*,D)  cannot  be  computed  until  R(*,D)  has  been 
computed. _  The  function  Rq(*,D)  represents  the  expected  damage 
when  the  defender  decides  not  to  use  his  second  volley  at  an 
unintercepted  RV. 


In  particular 


Rn(l,D)  =  1-s 


i.e.,  the  probability  that  the  target  is  destroyed  by  the  first 
RV  when  it  is  not  engaged.  We  now  show  that  Rn(*,D)  does  not 
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enter  in  the  determination 


D)  as  long  as  r<l,  and  there 


Q  f  Q  (  • 

need  not  be  computed.  The  result  is  intuitive:  If  you  have  a 
least  one  interceptor  left  to  engage  an  approaching  BY,  it  is 
better  to  use  it  immediately  (ever,  if  it  has  low  reliability) 
instead  of  saving  it  for  a  possible  subsequent  BY  (even  if  its 
first-shot  reliability  might  be  hi-gh). 

Lemma.  Assume  C:i_r<  1 .  Then  for  any  T>_1 


Proof .  Fix  D_>1.  We  will  use  induction  on  A. 

Since  Q(G,D)  =  0,  we  have 

{^(1,0)  =  rk(i-s)  for  any  k,  and 
so 

R(1,D)  <  1-s 

~  v  ■  A 

Qk(l,D)  =  qk  R( l,D-k)  <  1-s 

and  it  follows  that  (7)  is  true  for  A=l. 

For  any  k, 

Qk(A,D)  =  qk  R( A, D-k )  +  ( l-qk )  Q(A-l,D-k) 

=  qk[rk(l-s)  +  ( l-rk ( 1-s  )  )  Q(A-l,D-k)] 
+  ( l-qk  )  Q( A-l ,D-k ) 

=  qkrk(l-s)  +  ( l-qkrk( 1-s ) )  QCA-ljD-kO 
<_qkrk(l-s)  +  ( l-qkrk  ( 1-s  ) )  (A-l)(s-l) 


where  induction  is  used  to  get  the  last  inequality.  Thus 

Qk(A,D)  £  [q^r^  +  (l-q^r^C 1-s ) (A-l)] (s-1) 

and  the  quantity  in  square  brackets  is  clearly  bounded  above  by 
the  integer  A.  Thus 


Qk(A,D) 


<_  1— s  for  any  k 


and  the  result  follows 


Theorem.  If  0<_r,s<l,  then  RQ(*,D)  need  not  be  used  to  complete 
R  (  • ,  D ) . 

p roo f .  To  compute  R(*,D),  we  first  determine  n(D)  and  k**(D) 

from  ( 2 ) . 


;or  k=0,  we  have 


Rn ( A , D )  =  ( 1— s )  +  sQ(A-l,D) 


In  oarticular 


Rn(l,D)  =  1-s 


and  we  now  show  that 


Rn(l,D)  Rn ( A, D ) 


for  A=2 , 3 , . . . 


i.e.  , 


i.e.  , 


A(l-s)  >  (1-s)  +  s  Q(A-1,D) 


( A-l )  >  s ( A-l  +  Q ( A-l , D ) ) . 
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But  the  lemma  implies 


Q ( A-l , D )  i  ( 1-s ) ( A-l ) 


i  .  e  .  , 


(A-l)  >  s ( A-l  +  Q( A-l j D ) ) 


so  that  (9)  and  hence  (8)  is  true.  Thus 


max  R„(A,D) 

AeX+  i  1-s  =  R0(1.D) 


We  will  now  show  that 


R1(A,D) 


<  1-s 


so  that  k**(D)  of  (2)  is  not  0.  We  need  to  show  that 


r(l-s)  +  (l-r(l-s))  Q ( A-l , D-l )  <  A(l-s)  for  all  A>1 


but,  again  by  the  lemma, 


Q(A-1  ,D-1)  <  (A-l) ( 1-s ) 


so  that 


r(l-s)  +  (l-r(l-s))  Q( A-l , D-l )  <  r(l-s)  +  ( 1-r ( 1-s ) ) ( A-l ) ( 1-s ) 


«  [r  +  (l-r(l-s) ) (A-l)] (1-s) 


and  since  r<l,  the  quantity  in  square  brackets  is  strictly 
smaller  than  A  and  the  theorem  is  proven. 


Ex amp le .  With  (q,r,s)  =  (.1,.2,.3).  Table  1  exhibits  the  values 
Q(A,D)  for  A=0,l,...,10  and  D=0 , 1 , . . . , 10  .  The  first  column  of 
the  table  (Q(0,A))  represents  the  probabilities  of  destruction 
in  the  non-defended  case. 


Table  1.  Case  (q,r,s)  =  (.1,.2,.3) 


DESTRUCTION  PROBABILITIES  Q(A,D) 


A 

Q(A,0) 

Q(A,1) 

Q(A,2) 

Q(A,3) 

Q(A,4)  Q(A,5) 

Q(A,6) 

Q(A,7) 

Q(A,8) 

Q(A,9) 

Q(A 

t,io: 

0 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0. 

000 

1 

0.700 

0.070 

0.014 

0.014 

0.014 

0.014 

0.014 

0.003 

0.003 

0.003 

0. 

003 

2 

0.910 

0.721 

0.137 

0.033 

0.028 

0.028 

0.028 

0.017 

0.007 

0.007 

0. 

006 

3 

0.973 

0.916 

0.741 

0.199 

0.055 

0.042 

0.041 

0.031 

0.021 

0.012 

0. 

010 

4 

0.992 

0.975 

0.922 

0.760 

0.257 

0.081 

0.056 

0.045 

0.034 

0.025 

0. 

016 

5 

0.998 

0.992 

0.977 

0.928 

0.777 

0.311 

0.109 

0.079 

0.052 

0.039 

0. 

030 

6 

0.999 

0.998 

0.993 

0.978 

0.933 

0.793 

0.361 

0.176 

0.104 

0.060 

0. 

046 

n 

\ 

1.000 

0.999 

0.998 

0.994 

0.980 

0.938 

0.808 

0.418 

0.239 

0.131 

0. 

074 

3 

1.000 

1.000 

0.999 

0.998 

0.994 

0.981 

0.942 

0.825 

0.470 

0.296 

0. 

162 

9 

1.000 

1.000 

1.000 

0.999 

0.998 

0.994 

0.983 

0.948 

0.841 

0.518 

0. 

340 

10 

1.000 

1.000 

1.000 

1.000 

0.999 

0.998 

0.995 

0.984 

0.952 

0.855 

c. 

561 

THE 

NUMBER  TO  SHOOT  AT  NEXT  RV  IN 

1ST 

VOLLEY 

: 

0 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

IF 

MISS, 

THE  NUMBER  TO  SHOOT  IN  2ND  VOLLEY: 

0 

0 

1 

1 

1 

1 

1 

2 

2 

2 

2 

SLOPES (D) 

0.700 

0.361 

0.247 

0.190 

0.156 

0.134 

0.118 

0.105 

0.095 

0.086 

0. 

056 

Also  Included  In  Table  1  Is  the  firing  doctrine  which  yields 
the  Q(A,D)  values.  To  Incorporate  this  doctrine,  suppose,  for 
example,  that  the  defender  has  10  interceptors  at  the  target. 

If  he  sees  an  attacker  coming  in,  he  uses  one  of  these  10  in  the 


first  volley  and,  if  he  misses,  he  shoots  two  in  the  second  volley 
These  values  are  read  below  column  Q(A,10)  of  the  Q(A,D)  table. 

If  the  first  volley  is  unsuccessful,  but  the  second  succeeds,  the 
defender  now  has  10-1-2  =  7  interceptors  left  to  deal  with  any 
subsequent  attackers.  If  indeed,  a  second  RV  attacks,  the 
defender  again  uses  a  1-2  firing  doctrine  as  Indicated  below 
column  D=7.  If,  again,  his  first  volley  fails  but  his  second 
succeeds,  he  has  7-1-2  =  4  interceptors  left. 

A  possible  (but  unlikely*)  battle  history  with  D=10  is: 


RV  # 


1st  2nd  #  Interceptors 

Volley  Volley  Left  if  1st  Volley  Misses 


The  slopes  m(D)  of  the  lowest  maximum  damage  curves  which 
yield  an  upper  bound  on  the  expected  damage  are  also  given  in 
Table  1.  Thus,  for  example,  with  7  interceptors,  the  expected 
damage  Is  bounded  above  by  the  linear  function  F(A)  =  0.105A. 

Figure  2  exhibits  the  expected  damage  curves  for  increasing 
values  of  D. 

Often  one  wishes  to  determine  the  minimum  number  of  inter¬ 
ceptors  required  to  enforce  a  maximum  damage  function  of  a  given 
slope.  Figure  3  exhibits  the  curves  both  with  and  without  an 
SLS  capability  when  the  interceptor  failure  probabilities  are 
0.3  and  the  RV  failure  probability  is  0.1.  For  example,  if  one 
wishes  to  design  a  defense  whose  expected  damage  per  attacking 


*The  probability  that  this  would  occur  is 

(.lKl-.22)(.l)(l-.22)(.l)(l-.2)(.l)(l-.2)  =  5-9  x  10-5 


MTV 


Figure  2.  Expected  Damage  Curves  and  Maximum  Damage  Bounds  for  Varying 


weapon  is  0.15,  one  needs  6  interceptors  with  an  SLS  capability, 
and  ?  interceptors  without. 

Extensions  and  Implementation 

The  case  q=i  corresponds  to  no  SLS  capability.  The  case 
s=l  is  a  trivial  case  wherein  the  RV's  are  completely  unreliable 
sc  no  defense  is  needed. 

The  results  of  this  note  are  easily  extended  to  allow  for 
additional  volleys. 

A  FORTRAN  program  has  been  coded  to  generate  the  values 
5(A,T)  for  any  range  of  A,D  values. 


*  .  a 


Here  we  address  the  case  where  the  attack  is  "sequential" 
i.e.,  there  is  enough  time  between  successive  attackers  that 
can  be  ordered  and  the  attack  size  is  known.  We  define,  as 


oe: o re : 


?(A,D)  *  probability  that  the  target  survives,  given 
that  it  is  under  attack  by  A  missiles  and  is 
optimally  defended  by  D  defenders. 


Obviously,  if  the  defender  knows  the  value  of  A,  he  will  defend 
uniformly  according  to  the  result  of  Appendix  A.  (A  simultaneous 
attack  can  be  considered  sequential  by  numbering  the  attackers  in 
any  order.  ) 

However,  if  the  defender  has  a  shoot-look-shoot  capability, 
and  sufficient  time  between  arrivals,  he  can  choose  to  structure 
his  defense  in  volleys,  with  the  prospect  of  saving  defenders  for 
use  against  future  attackers. 

Suppose  the  defender  has  time  for  two  volleys  against  each 
incoming  attacker.  Let  a  be,  as  before,  the  kill  probability  of 
an  attacking  missile.  Let 


d  =  probability  that  a  defending  Interceptor  will  destroy 
an  attacking  missile  in  the  first  volley 


e  *  probability  that  a  defending  interceptor  will  destroy 
an  attacking  missile  in  the  second  volley. 


d(A)  *  number  of  interceptors  to  shoot  at  the  first  of  A 
attacking  missiles  in  the  first  volley 


1  Source:  Appendix  C  of  Reference  ( 1  ]  on  page  R- 1 . 


e(A)  =  number  of  interceptors  to  shoot  at  the  first  of  A 

attacking  missiles  in  the  second  volley,  given  that 
the  first  volley  has  failed. 

Then 

l-d(l-d)d(A)  is  the  probability  that  the  first  volley  is 
successful, 

( 1-d ) d ^ A ^ ( 1- ( 1-e ) e ^ A ^ )  is  the  probability  that  the  first 

volley  fails  but  the  second  is 
successful , 


( 1-d ) d^ A ^ ( 1-e ) 6 ^ A ^ ( 1-a )  is  the  probability  that  both 

volleys  fail  and  the  attack  also 
fails. 

The  following  recursion  holds: 

P ( A, D )  =  max  | (l-(l-d)d(A) )  P(A-l,D-d(A) ) 


d  ( A ) ,  e  ( A )  e  I 
d  ( A )  +  e  ( A )  <  D 


+  (l-d)d(A)(l-a(l-e)e(A))P(A-l,D-d(A)-e ( 


with 


P ( 0 , D )  *  1  for  all  D  e  I+ 


Given  a,  d  and  e,  the  recursion  can  be  solved  by  dynamic 
programming  to  determine  P(A,D)  =  p^ .  Note  that  the  solution 
of  this  recursion  would  agree  with  the  results  of  Appendix  A 
in  the  case  where  e*0. 


Obviously,  the  above  recursion  could  be  extended  if  the 
defender  had  more  than  two  opportunities  to  protect  himself. 


DISTRIBUTION 


IDA  PAPER  P-1875 

ROBUST  PREALLOCATED  PREFERENTIAL  DEFENSE  MODEL 


Copies 


Director, 

Office  of  the  Joint  Chiefs  of  Staff 
Washington,  D.C.  20301-5000 
ATTN:  Director,  Joint  Analysis  Directorate 

Defense  Technical  Information  Center 
Cameron  Station 
Alexandria,  Virginia 

Deputy  Chief  of  Staff  for  Research  and  Development 
and  Acquisition 
Department  of  the  Army 
Room  3A474,  The  Pentagon 
Washington,  D.C.  20310 
ATTN:  Director,  Missile  and  Air  Systems  Division 
RDA/DAMA-WS 

Office  of  the  Chief  of  Staff 

Department  of  the  Army 

Ballistic  Missile  Defense  Program  Office 

P.O.  Box  15280 

Arlington,  Virginia  22215 

ATTN:  DACS-BMZ 

Commander 

Department  of  the  Army 

Ballistic  Missile  Defense  Systems  Command 

P.O.  Box  1500 

Huntsville,  Alabama  35807 

ATTN:  Library 

Headquarters 

Department  of  the  Air  Force 
Assistant  Chief  of  Staff  Studies  and  Analysis 
Room  1E388,  The  Pentagon 
Washington,  D.C.  20330 
ATTN:  Library 

Directorate  of  Aerospace  Studies 
Deputy  Chief  of  Staff,  Plans  and  Programs 
Headquarters,  Air  Force  Systems  Command 
K inland  AFB.NM89117 
ATTN:  Library 


Copies 


The  Rand  Corporation 
P.O.  Box  2138 

Santa  Monica,  California  90406-2138 
ATTN:  Library 

The  Rand  Corporation 
2100  M  Street,  N.W., 

Washington,  D.C.  20037 
ATTN:  Library 

Hudson  Institute,  Inc. 

Center  for  Naval  Analysis  (CNA) 

P.O.  Box  11280 
Alexandria,  Virginia  22311 
ATTN:  Library 

Los  Alamos  National  Laboratory 
P.O.  Box  1663,  Mail  Station  5000 
Los  Alamos,  NW  87545 
ATTN:  Library 

University  of  California 

Lawrence  Livermore  National  Laboratory 

P.O.  Box  808 

Livermore,  CA.  94550 

ATTN:  Library 

Analytic  Services,  Inc.  (ANSER) 

Crystal  Gateway  3 
1215  Jefferson  Davis  Highway 
Arlington,  Virginia  22202 
ATTN:  Library 

Teledyne-Brown  Engineering 
Cummings  Research  Park 
Huntsville,  Alabama  35807 
ATTN:  Library 

Teledyne-Brown  Engineering 
1250  Academy  Park  Loop 
Colorado  Springs,  Colorado  80910 
ATTN:  Library 

SAIC 

1710  Goodridge  Drive 
McLean,  Virginia  22102 
ATTN:  Library 

McDonnell-Douglas  Astronautics  Company 
5301  Bolsa  Avenue 
Huntington  Beach,  CA.  92647 
ATTN:  Library 


vi 

Vi 


Sparta,  Inc. 

4901  Corporate  Drive,  Suite  102 
Huntsville,  Alabama  35805 
ATTN:  Library 

System  Planning  Corporation 
1500  Wilson  Boulevard 
Arlington,  Virginia  22209 
ATTN:  Library 

Institute  for  Defense  Analyses 
1801  North  Beauregard  STreet 
Alexandria,  Virginia  22311 


N  , 

ATTN: 

s 

Dr.  J.  Bracken 

2 

Dr.  W.J.  Schultis 

1 

Mr.  R.B.  Pirie 

1 

Mr.  S.J.  Deitchman 

1 

•  * 

Control  &  Distribution  5 

i 


! 

>1 


2 

DL-3 

I 


V.V.V.'V'.V 


